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Abstract 

The first flights of the NASA/Marshall airborne C02 Doppler lidar wind 
measuring system were made during the summer of 1981. Successful 
measurements of two-dimensional flow fields were made to ranges of 15 km 
from the aircraft track. This report examines the characteristics of the 
data obtained, and summarizes a study of various artifacts introduced 
into the data set by incomplete compensation for aircraft dynamics. Most 
of these artifacts can be corrected by post processing, which reduces 
velocity errors in the reconstructed flow fields to remarkably low 
levels. 
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As part of the NASA severe storms program NASA/Marshall has * 

constructed an airborne lidar system designed to measure two-dimensional 
horizontal flow fields in the atmosphere. This system uses a pulsed 
coherent C02 laser operated much as a DoppLer radar. Velocity components 
along the lidar beam are measured by observing the Doppler shifts in the 
signals returned from naturally-occurring aerosol scatterers. The beam is 
steerable within a forty-degree cone to the left of the aircraft; by 
taking measurements with the beam oriented in the two horizontal extremes 
of the cone (70 and 110 degrees to the left of the aircraft heading) two 
components of the horizontal flow field are sensed* From these two 
components it is possible to reconstruct the two orthogonal components of 
horizontal flow. 

This complex system was assembled by a group of NASA personnel and 
contractors over a three-year period. The first data flights using the 
NASA/Ames CV990 aircraft were completed during the summer of 1981. Data 
was obtained in a variety of situations: boundary-layer and orographic 
flows were observed in California, Colorado, Oklahoma, Nevada and. 

Montana, while flows in the vicinity of convective storms were observed 
at many levels in Montana during participation in the CCOPE experiment. 

It is no simple matter to operate a Doppler radar using an aircraft 
as a platform. The velocity of the aircraft is high compared with the 
atmospheric velocities which must be measured. This motion of the 
measuring instrument must of course be removed from the measurement, a 
feat which can only be accomplished by very accurate knowledge of the 
platform velocity and all the relevant angles relating aircraft motion, 
attitude and instrument pointing. 

This much is required for accurate measurement of the velocity 
component parallel to the lidar beam. To obtain horizontal flow vectors, 4 
a second component must be measured by steering the lidar beam as 
described above. During a measurement run a series of measurements is 
obtained at each observation angle. Each series consists of a grid of 
velocity component measurements typically separated in range and track by 
300 ra (the range gate width and the scanner period). The two grids of 
measurements are then superimposed and processed to extract the velocity A 
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vectors, just as vectors are derived in ground-based multiple Doppldr 
measurements . 

To make such interpretation possible the relative registration of 
the measurement points must be sufficiently well known - within a 
fraction of a grid spacing - to permit correlating the proper 
observations with each other, a requirement which places stringent 
limitations on the allowable errors in aircraft navigation and lidar beam 
pointing. 

The flow-vector estimate produced by operations on the two component 
measurements can, at best, be only as good as the original radial 
velocity measurements. In fact, it is worse, since the poor geometry of 
the twp measurements (with only a 40-deg included angle) triples the 
error in the vector component parallel to the aircraft track; such 
geometry mandates both very high accuracy on the part of the Doppler 
estimator and very high stability on the part of the lidar local 
oscillator and transmitter. 

The state-of-the art in aircraft navigation and lidar systems is 
strained by such stringent accuracy requirements. Therefore the success 
of the 1981 field tests is all the. more remarkable, and is a tribute to 
the efforts of all concerned. 

While the program was a success in the sense that useful 
measurements of several types of flow fields were obtained, many of the 
measurements are clearly contaminated. The purpose of this report is to 
survey the problems present, to suggest corrections where corrections are 
possible, and to describe problems not now understood. 

Most of the artifacts present in the data sets appear to stem from 
incomplete compensation for aircraft dynamics. Since the motion of the 
platform is so critical to these problems, Section II is devoted to 
examples of time-series plots of several important aircraft parameters: 
pitch, roll, velocity, etc. Examples of radial-velocity measurements and 
lidar pointing are also presented. 

Section III uses "real-time" plots (obtainable on the aircraft in 
essentially real time for use in experiment management) to demonstrate 
the characteristic types of artifacts found in the raw data sets. 

Section IV is devoted primarily to the major cause of these 
artifacts - incorrect interpretation of the aircraft drift angle. This 


m 
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section opens with time-series examples of the important parameters and 
their derivatives. The cause of the drift-angle problem is suggested, and 
the argument is reinforced by means of simulated uniform flow fields 
wherein artificai drift-angle errors have been introduced. A first-order 
correction is suggested and tested - first on a time series of mean 
radial velocity, and then upon actual, flow fields. This correction 
removes about 75% of the largest source of error in the data. 

Other sources of error are present, however, as discussed in Section 
V. Included in this category are correctable errors due to incorrect 
horizontal pointing information, uncorrec table errors in the elevation 
angle of the lidar beam, and errors due to lidar moding and 
local-oscillator drift which are not yet understood. 

Section VI is devoted to examples of interesting flow fields 
observed during the 1981 program. These samples (corrected to first-order 
for drift-angle errors) demonstrate the potential of the system for 
investigating orographic flows, boundary layer flows during convection, 
gust-fronts in clear air, and clear-air flow in the vicinity of 
convective storms at mid- and upper-levels. 

The detective work in unravelling the mysterious problems with these 
data sets is not yet complete. It is felt that a significant part of the 
remaining errors in the data is correctable, and techniques for such 
correction have been devised but not tested. Section VII discusses 
briefly this work in progress. 

During investigation of this mass of field test data the author 
attained a certain degree of familiarity with it, and a good deal of 
software was written to allow examination of the data and testing of 
various correction algorithms. While this software was written strietly 
for internal use, it does constitute a complete set of software for 
processing the data from original tape to gridded flow fields. It was 
thought worthwhile to include the software in this report, not so that 
all of the software will be used as such, but rather to explain, in a 
definitive way, questions of treatment. 

Appendix A lists the format of the raw data tapes, with some 
comments on problems in the data headers. Appendix B lists the software 
just mentioned, with enough description to orient the reader. 
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As indicated in the introduction, aircraft dynamics are critical 
factors in the error budget of the data sets. The data acquisition system 
attempts to monitor the relevant aircraft parameters and to correct for 
platform motion and attitude, Correction is accomplished in two ways; the 
lidar scanner attempts to compensate for attitude changes, and the lidar 
second local oscillator attempts to compensate for velocity changes. 

Section III will discuss the degree to which these corrections are 
successful, Th« intent of this section is to provide examples typical of 
platform motion during the field tests. The data set selected is run 10 
of flight 19, a run at low altitude up the western side of the San 
Joaquin Valley in California. This run was selected because of its 
unusual length (about 40 minutes) and because the rather uniform wind 
field allows measurement errors to be easily seen. 

Figure 1 shows ground speed as estimated by the inertial navigation 
system (INS) as a function of time over the 40-minute data collection 
period. The resolution of this measurement is 1 kt, and it is subject to 
errors due to INS drift. Since a fraction of this velocity of about 
cos(70 deg) must be removed from the measurement to correct for platform 
motion, it is clear that the ultimate accuracy of an individual Doppler 
measurement of radial velocity referenced to ground cannot exceed 0,17 
m/s = 1 kt * cos (70). 

The variation of altitude during the run is shown in figure 2. The 
aircraft normally flies at constant pressure altitude, to an accuracy of 
about 15 meters. Complications arise for the flow field measurement 
whenever the altitude is changed; at the longer ranges the delay between 
forward and aft measurements at the same point may be up to one minute, 
and if the altitude has changed the two measurements will reflect flows 
at two separate altitudes. This is significant in the boundary layer. 

Figure 3 shows aircraft true heading during the run. The general 
trend of the plot reflects the aircraft course along the curving western 
edge of the San Joaquin Valley. The fine structure, amounting to about 
one degree peak-to-peak, reflects aircraft control-system stability in 
mild turbulence. The oscillations of a few degrees are aircraft responses 
to changing crosswind. True heading is obtained from the INS navigation 
unit with a resolution of 0.44 deg. 
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True airspeed is used by Che data acquisition system only to 
estimate the probable Doppler shifts in the data (since the nominal 
Doppler shift is expected to be TAS*cos(7Q) ). The form shown in figure 4 
is quite similar to ground speed (figure l) , with a small offset due to a 
headwind. 

The pitch of the CV990 aircraft is rather stable, as shown in figure 
5, with transients on the order of 1 deg whenever altitude or airspeed ia 
altered. Pitch is obtained from the INS unit with a resolution of 0.44 
deg., and is used by the lidar scanner to hold the lidar beam in the 
horizontal plane. Aircraft roll, shown in figure 6, 1 is used by the 
scanner in the same manner. Roll is subject to excursions on the order of 
degrees, often with a well-defined period. 

Aircraft drift angle is critical to the correction for platform 
velocity. The plot shown in figure 7 exhibits a typical short-term 
variability of about 1 deg peak-to-peak, with larger variations where the 
flow field is complex. Drift angle often Shows resonant oscillations 
related to the feedback characteristics of the aircraft control loops. 
One ; ogsin the resolution of this measurement is 0.44 deg. 

Using knowledge of ground speed, airspeed, and drift angle the INS 
unit estimates the wind vector in the vicinity of the aircrafts Figures 8 
and 9 show the magnitude and direction of this estimate for the run. 
While this measurement is not always reliable, it has been a useful 
comparison for the lidar measurements. 

Examples of the raw lidar radial velocity measurements are shown in 
figures 10 and 11. The radial velocities have been averaged over 30 range 
gates. While these measurements cannot be compared directly with the INS 
wind estimates (since the INS estimate is a magnitude, and the lidar 
measurements are components), several features on the curves agree and 
demonstrate that the lidar system is measuring something related to the 
wind field. Note that the fine structure in the lidar measurements is on 
the order of 1 m/s . 

The look-angles at which the forward and aft lidar measurements were 
taken are shown in figures 12 and 13. The form of the curves closely 
approximates the aircraft true heading (figure 3), with 70- and 110-deg 
offsets. The resolution of these angular measurements is 0.1 deg. 
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Figure 5: Pitch angle 
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Figure 10 : Forward radial velocity 
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Figure n : Aft radial velocity 




Figure 12 ; Forward beam pointing angle 






Figure 13: Aft beam pointing angle 




It should be noted that the examples, shown above are typical of 

* tt 

flight in mild turbulence. At higfct attitudes aircraft parameters are 
sometimes maintained to tighter tolerances, while in more severe 
turbulence excursions can be markedly higher. 


III. Raw data characteristic# 


The data acquisition so£tware used during the 1981 flight program 
provided £or ’’real-time" plots of horizontal flow fields. The plots were 
updated on a scan-by-sean basis to allow rapid feedback to the 
experimenters. This section uses six of these plots to illustrate the 
characteristics of the uncorrected data. The cases were chosen for their 
relatively uniform flow fields, which make data errors more visible. 

Figure 14 illustrates a rather smooth flow field with few artifacts 
present. As such it is typical of the better measurements . The format is 
similar to that of the figures to follow: the flight track is at the 
bottom of the figure, left to right. The grid spacing of the measurements 
is roughly 300x300 m, with a nearest^neighbor criterion used to match 
forward and aft scan points. At longer ranges, corresponding to the top 
of the figure, useful returns were not obtained, as indicated by various 
letter codes. 

Figure 15 is taken from earlier data on the same flight. The flow 
vectors are confused and difficult to interpret. One might infer that 
there is a complex flow pattern arising from convection present; in fact, 
errors to the extent of a few m/s are present in the vectors, and the 
pattern of the errors is not readily apparent. Near the top of the figure 
are a few "wild" measurements, often found at long ranges where the 
returns are weak. The apparent complexity of this measurement on a nearly 
uniform field should be a warning to interpreters of flow fields: make 
sure the data is significant before engaging In deep interpretation. 

More obvious artifacts are present in the example shown in figure 
16. There is a definite tendency for vectors to line up on the 70- and 
110-deg radials at which lidar measurements were made. This triangular 
appearance of the data was one of the first characteristics noted on many 
of the early plots. Also seen in this figure, in addition to a few 
isolated "wild" measurements, is a row of uniform, strong vectors in th/e 
upper left. A rather frequent a result of multiple-mode interference in 
the laser system, they are easily rejected. 

This same triangular effect takes on a different appearance when the 
flow is primarily parallel to the aircraft track, as in figure 17. Note 
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that there is a modulation of the vector length in a roughly triangular 
pattern. Once again there is a line of strong, constant vectors due to 
laser moding. 

When the system errors are periodic the errors seen in the previous 
examples take the form of braiding of the flow vectors, as in figure 18. 
One might imagine that the flow is dominated by organized horizontal 
rolls, or perhaps influenced by canyon topography; but the effect is 
instead produced by aircraft dynamics, as discussed in Section IV. 
Another example of the same effect is shown in figure 19, where the data 
was taken at 12,400 feet and the structure is due almost entirely to 
aircraft dynamics. 

The moral of this section is clear. Some very interesting effects 
appear in the data, but unfortunately they are instrumental in origin. In 
many cases the artifacts are so severe that they completely obscure the 
nature of the flow field. This is particularly true in the boundary 
layer, where the flow-field perturbations of interest in convection 
studies are quite small. 
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This section is devoted to an explanation of the most significant 
source of artifacts in the data* erroneous drift-angle information. To 
establish conventions for various angles, refer to figure 20. The body of 
the aircraft is aligned along the vector THDG (true heading), but the 
center-of-gravity of the aircraft is travelling along the vector TRK 
(track) relative to the ground. The two vectors do not coincide due to 
the assumed presence of a wind from the left; the difference between the 
two vectors is the drift angle, which in this case is positive. The 
direction of the lidar measurement is given by vector LOS 
(line-of-sight), and this angle is measured by angle SCAN with reference 
to THDG. 

Note that the ground velocity is along vector TRK, It is the 
component of this velocity parallel to LOS which must be subtracted from 
the data measurements in order to reference the wind flow to the ground. 
Thus the angle between the laser measurement and the aircraft motion is 
SCAN+DRIFT . The angle SCAN is closely controlled and monitored by the 
lidar scanner; thus the angle DRIFT is crucial for proper application of 
the velocity correction. 


WIND 



DRIFT 


Figure 20: Lidar-aircraf t geometry 
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With this in mind one can conclude that drift angle and aircraft 
ground speed are the critical parameters required for removal of the 
component due to aircraft motion. Three aspects of these two parameters 
can cause errors in this correction: resolution, accuracy, and time 
delay. Clearly if the resolution of the measurement is limited the 
correction will be correspondingly limited. Similarly, if .there are 
errors in the measurement of drift angle or ground speed there will be 
errors in the correction. Finally, if the correction quantity is out- 
dated when it is used there will be errors, provided that the quantity 
changes fast enough. The latter effect dominates errors in the data sets. 

Figure 21 shows a 2-minute sample of raw radial-velocity 
measurements taken during a portion of run 10, flight 19. The two plots 
are for the forward and aft measurements, and 30 range gates have been 
averaged together. The wind field was relatively uniform. The most 
obvious features of the figure are excursions on (She order of a few m/ s 
which are correlated in the two plots. 0n\'? can easily show that such 
latge excursions in a mean velocity are unlikely: the velocities are 
averaged over 10-km in range, and show changes in tha means of up to 4 
m/s in 300-m of flight. This implies a shear of 0.01/s over 10 km, an 
order of magnitude higher than that typical of the turbulent boundary 
layer, in a situation which was relatively smooth. Further, rapid changes 
in the two measurements are not consistent with correlation between them, 
since the averaged regions diverge from the aircraft, and are separated 
by 7 tggg at 10-km range. We may safely conclude that most of the 
excursions in figure 21 are due to some form of measurement error. 

A study of the continuity in range of individual lidar shots 
indicates that the Doppler estimation errors within each shot are quite 
low. The source of error is thus likely to lie in the correction of the 
measured velocities to ground-based coordinates. Therefore the next step 
is to investigate these suspect measurements. Figure 22 shows two such 
measurements for the same 2-min time period. The relationship between 
true heading and the "errors" of figure 21 is not apparent, but clearly 
drift angle is somehow correlated. 
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Figures 2? and 24 plot the derivatives of true heading and the two 
iine-of-sight angles of the forward and aft scans. These derivatives, 
which are estimated by taking 2-scan (approximately 2.2-sec) differences, 
ar>; also somewhat correlated with the "errors" in figure 21. 

Finally, figure 25 shows the 2-scan derivative of drift angle, 
plotted separately for the forward and aft scans. If this figure is held 
face to face with figure 21 a remarkable correlation will be observed. 
Not only is the correlation very accurate, but there is a time delay 
between the two figures amounting to about 2 seconds. 

Obviously, since there is something remarkably similar between 
drift-angle changes and apparently erroneous mean radial velocity 
measurements, the next step was to find the connection. Simulation was 
used to ascertain the impact of oscillatory drift-angle errors upon 
derived wind fields, given a variation in radial velocity of a few m/s 
and differences in drift angle on the order of 1 deg. The results are 
shown in figures 26-28. In each case a uniform wind of 10 m/s was 
assumed, and an artificial sinusoidal perturbation in drift angle of 
about 0.3 deg was used in the model. 

In the case of figure 26 the assumed constant wind field was normal 
to the aircraft track (along the bottom of the page in these three 
figures). The effect of the drifjt-angle perturbation is substantial, even 
though the perturbation was only 0.3 deg. When the wind field is parallel 
to the aircraft track as in figure 27 the appearance is quite different: 
a modulation of the vector length very much like that seen in figure 17. 
Finally, with the wind vector at 45 deg relative to the aircraft track 
the effect is a braiding of the flow vectors identical to that seen in 
figures 18 and 19. 

To quantify this effect note that in terms of the geometry of figure 
20 the magnitude of the velocity correction is GS*COS(SCAN+DRIFT) . 
Assuming that the quantities GS (ground speed) and DRIFT are subject to 
errors (from whatever source) EGS and EDR, one can expand the expression 
and determine the errors in the correction term. The expansion contains 
three terms: 

EDR * GS * SIN ( SCAN+DRIFT ) 

EGS * SIN (SCAN+DRIFT) 

EDR * EGS * COS (SCAN+DRIFT) 
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Figure 23: 2-lag differential true heading 
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Figure 25 Forward and aft differential drift angle 





Figure 28: Simulated 15-deR flow 
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Typical value? for GS and SCAN+DRIFT are 13 0 m /$ and 70 deg. In the plots 
shown above typical values for EGS and EDR are 1 m/s and 2 deg, for 2-sec 
differences; these values yield estimated contributions for the three 
terms of 4, 0.3 and 0.03 m/s, respectively. Clearly the cross term is not 
significant, and the second term is relatively small. However, the first, 
term is quite significant, and produces errors of the proper magnitude. 

With this confirmation of the critical factor a correction of the 
error can be attempted. Using the 2-lag drift difference as the error, 
and the geometry just described, the mean radial velocities can be 
corrected as shown in figure 29. The data is the same as shown in figure 
21, but the variance has been reduced substantially. 

Before further discussion of correction algorithms an explanation of 
the source of the drift-angle error is in order. The error arises not 
because the drift angle is in error, but because the drift angle used for 
velocity subtraction by the data acquisition system is not the current 
drift angle. The data gathering process is a complex one, depending upon 
the computer in the inertial navigation system, the computer in ADDAS 
(the CV990 data gathering system), and the computer in the lidar system. 
To take the case of drift angle -in particular, it is first calculated by 
the INS system in what is termed the "slow" loop, updated at 0.9-sec 
intervals. This data is then output on the BCD bus at 0.9-1.0-sec 
intervals adding (on the average) an additional 0.48 sec to the delay. 
When this data is read by the lidar computer at the start of scanner 
motion it is about 1.38-sec old. By the time the scanner has moved to 
target position and the midpoint of the data sample has been reached 
another 0.84 sec has elapsed, for a total delay of 2.22 sec. 

This 2.22-sec delay is nearly equal to 2 scan periods (2.2 sec on 
the average for this data set), and this explains why the 2-lag 
drift-angle difference correlates so well with the mean radial velocity 
"error". The 2. 22-sec figure is a theoretical one, and must be confirmed 
for each situation. A routine is included in program FILED1 (see Appendix 
B) for finding the optimum 4~point transversal filter (or. interpolator) 
for this drift-angle correction. For most of flight 19 the optimum filter 
coefficients were found to be 0.32, 0.56, 0.12 and -1.0, corresponding to 
a time lag of 2.42 sec. 
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Using this optimal filter the data oE Eigura 21 has been corrected 
to produce figure 30. The results are reasonably smooth) and definitely 
superior to the simpLe 2-lag result of figure 29. The wind field produced 
by the uncorrected data of figure 21 is shown in figure 31, and the 
corrected version is shown in figure 32. Considerable improvement is 
evident. 

As an additional example another raw radial velocity data set is 
shown in figure 33. The velocity excursions are marked and regular) and 
appear to be due to an oscillation in the aircraft control system. 
Applying the optimal drift-angle correction to this data set produces 
figure 34. The comparable raw and corrected flow fields for these 
examples are shown in figures 35 and 36. The braided structure seen 
earlier is present in the former figure) but nearly gone in the latter. 

As a final example a situation which exaggerates drift-angle effects 
is shown in figure 37. The aircraft began a slight turn midway through 
the example) producing artifacts consisting of counter-rotating flows in 
this uncorrected plot. 

While drift-angle delays have been identified as the major 
contributor to the artifacts evident in the data sets, other effects are 
present as well, and the algorithm for correction discussed above is not 
completely satisfactory even for drift-angle correction. These additional 
considerations are discussed in the next section. 
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Figure 33; Uncorrected mean forward and aft lidar radial velocities 




Corrected mean forward and aft lidar radial velocities 
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Figure 35: Uncorrected 
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Figure 36: Corrected flow field 
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Figure 37: Flow field with errors due to aircraft turn 





V. Other data errors 


V.1 Errors due to delays in INS data 

i 

As mentioned in Section IV the correction for platform motion is 
also sensitive to delays in the INS ground-speed estimate. The error from 
this source is on the order of 0.3 m/s; though smaller than that due to 
the delay in drift angle it is still significant and should be taken into 
account • 

The scanner control computer also relies upon INS attitude 
information in calculating the beam-pointing corrections required to 
compensate for roll and pitch. Errors due to the delayed angular 
information take the form of horizontal and vertical pointing errors on 
the order of 1 deg. The horizontal errors are easily corrected by 
interpolation of the INS estimates of true heading) with slight 
corrections from roll and pitch interpolations. Correction for the 
vertical errors is not possible. The magnitude of the errors may be 
estimated by interpolating the roll and pitch readouts) but the fact that 
the lidar beam did not travel in a horizontal plane cannot be changed. A 
1-deg error in the vertical translates into a 175-m vertical displacement 
at 10-km range; this distance can be quite significant in the boundary 
layer. An example of this type of error (not included in the figures) is 
a case where the aircraft was flying about 100 m above the boundary layer 
in very smooth air; the measured wind field was very smoot/i with the 

exception of a few shots where the lidar beam dipped into the boundary 
layer. Since shear at the boundary-layer interface was strong the 
resulting flow field was discontinuous. 

The problem of estimating and correcting for these data-delay errors 
is complicated by the fact that the lidar system does not operate at a 
constant duty cycle. The time between scans may vary substantially, 

depending upon the mode chosen by the operator. Thus the data set may 

contain a series of scans separated by 1.1 sec, with scans requiring over 
2 sec interspersed. Correction formulas used thus far have assumed 
constant scan rates, and more complex algorithms must be developed to 

allow for the variable scan rate. 

t . \ 

V.2 Hardware problems 


Reference has bean made above to a problem in the Udar system 
termed M moding u » The laser is capable of oscillation in more than ona 
mode, and whan it doa» so tha modes mix on tha photodetactor and create 
Strong monochromatic signals* Such a case is evident in figure 17 as a 
row of strong, constant vectors. Data in such cases is lost, but it is 
not diffiee.lt to recognise this type of error, 

Many of the low-level flights in California were near foothill 
regions, In several cases Udar returns from terrain were obtained* 
allowing cheating of velocity correction (since ground returns should 
shown aero velocity). Tha actual location of the aircraft could be 
obtained vary accurately from pictures taken periodically during the 
flights by a downward-looking camera. When terrain returns were compared 
with topographic maps it waa discovered that there may be some fixed 
errors in the vertical elevations of the forward and aft beams. In 
particular, returns were obtained from the aft beam which indicated a 
bias of about -0,4 deg in that beam (that is, the tended fee he below 
the horizontal). No terrain returns were obtained from the forward beam 
in situations where terrain was 1,0 deg above the horizontal* indicating 
a bias of 1,0 deg or more for the forward beam, Terrain returns have also 
suggested tha possibility of a small bias error in Udar range* 

Other hardware problems are suspected but less easily seen. There 
are cases in the data where velocity measurements appear to be accurate 
near the aircraft but deteriorate with range, Such cases are sometimes 
due to errors in the vertical pointing of the beam, but it appears 
possible that some are due to frequency modulation of the 
local-oscillator laser, This laser is the master-oscillator laser for the 
transmitter as wall} as such, its frequency is by definition correct at 
the time of the pulse transmission* but if its frequency is not stable 
this relationship will deteriorate with time (or range). Strong 
mechanical vibrations associated with flight through turbulent air may 
excite such frequency drift. The extent to which this error source is 
significant is not known, 

Evident in the data are yet other types of errors whose sources are 
not known, fortunately most of these errors are large errors, and the 
editing software developed previously can affectively reduce their 
effects , 


o 


.-..a ^ .1..-..-:.-........... : a s^:daOSte^t... . 

VI, yKndflow example# 

this section presents examples of several types of Clows observed 
during the 1981 flight tests* These examples have been corrected to first 
order Cor drift-angle and true-heading delays* but they have not been 
smoothed . 

Figure 38 shows a very uniCorra Clow observed above the boundary 
layer at 9000 ft» The deep boundary layer was just below the aircraft* 
and it is likely that the discontinuous measurements at the right-center 
in the figure are due to the lidar beam penetrating the boundary layer. 
The discontinuous measurements at the left-center of the figure may be 
due to laser moding* since they persist down to aero range. 

Flow in the vicinity of a cumulus cloud is shown in figure 39. The 
regions in the figure with no measurements represent cloud regions not 
penetrated by the lidar. The velocity field is measured in the clear air 
around the cloud* and probably at the outer edge of the cloud itself. In 
some cases of this type regions at the sides of the cloud are not 

measured, since either the forward or the aft lidar returns were shadowed 

by the forward edges of the cloud, 

’ Numerous examples of orographic flow modification are prosent in 

data taken in California at low levels. Figure 40 shows south-easterly 
flow impinging upon the foothills of the Sierra Nevada near Fresno. As 
the air mass rises against rising terrain the lidar sees flow 
progressively nearer the ground, where the flow is modified by the local 
topography. Figures 41 and 42 show flow diverging from the Carquinez 
Strait into the Sacramento-San Joaquin Delta. In figure 41 the southern 

portion of the Clow exits from the strait at the lower right, forming 

eddies in the foothills in the upper part of the figure. In figure 42 the 
northern portion of the flow curls around a mountain forming the terminus 
of the strait. 

Measurements were made in the San Gorgonio Pass region of southern 
California. At this point strongly divergent flow exits the narrow pass, 
offering opportunities for the operation of wind turbines. Figure 43 
shows this strong flow entering at the lower right; weak confused flows 
are evident over foothills at the upper left. Figure 44 shows flow 
somewhat farther beyond the mouth of the pass, where complex patterns 
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Figure 40: Valley flow near Sierra Nevada 
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Figure 41: Southerly flow from Caruulnez Strait 
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Figure 43: Exit flow from San Gorgonio Pass 
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Figure 44: Exit flow from San Gorgonio Pass 






result from flow over terrain features (Including a small hill at lover 
center) * 

Some of the most spectacular flows observed were at lov levels 
beneath convective storms. Figure 45 shows at first sight an 
uninteresting flow field. The two disturbed areas are quite significant, 
however: they represent the horizontal signatures of areas of convection 
or outflow* A more dramatic example is shown in figure 46. A relatively 
uniform low-level flow is perturbed by an outflow descending from a 
convective structure at the upper Left. The structure itself is not seen, 
as the lidar returns are apparently attenuated by cloud in that area. 

Still stronger outflow - a true gust front - is shown in figure 47. 
The flow seen at the lower left is typical of a large area now shown in 
the figure, and may be taken as the unperturbed low-level flow. The 
strong flow crossing most of the figure is outflow from a convective 
structure of considerable size beyond the top of the figure. Similar 
merging of outflows with the low-level flow is shown in figures! 48 and 
49. Finally, in figure 50 a dramatic interaction between the outflow at 
the top of the figure and the low-level flow is apparent. * 
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Figure 46: Small gust front entering uniform flow 







Idg 251 — ► 1500' AGL 333-m Grid 10 m/. 


r 


4 > 

8 

L 

< 4 - 

4 > 

• 

0 

U) 


2 

0 

4 > 

c 

0 


00 

ca 

i 

N 



oHui;y 


6 


0 > 

u 

3 

u 

Pn 


- 68 - 


-f root interaction 







Montana qu®t front Hdq 078— ♦ 1300’ AGL 500-m Grid 10 m/® 



aHiiuy 


- 69 - 


Figure 48: Gust-front interaction 
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Figure 49: Gust-front interaction 
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Figure 50: Gust-front interaction 



VII* Continuing work 


1. Correction of navigation-delay errota. An optimal filter has been 
implemented for the drift-angle case* Similar filter# need to be 
developed for ground-speed correction, lina-of-sight angle r.ofrection in 
the horizontal, and line-of-sight angle estimation in the vertical. These 
filters must take into account the variable spacing of the scans in time. 

2. Continuity-argument filter. After all known sources of error in 
the mean radial velocity measurements have been estimated and corrected 
Cor, there will remain a degree of uncertainty in the moan radial 
velocities. This will be due to errors in the INS estimates of attitude 
and velocity, and to the random nature of the delay between INS 
measurements and INS outputs. These errors are inherently unknowable, but 
they can for the most part be corrected by an argument involving 
windfield continuity. 

Consider a time-series of moan radial-velocity measurements. Each 
item in the time series is an average over (for example) lO-km of range. 
The variance of each item will be due to that portion of the velocity 
fluctuation spectrum which Is not averaged out by the 10- km integration. 
That is, the major part of the variance of the true time series will be 
due to spatial variations in the wind field on the order of several km. 
The observed time series will contain additional variance due to the 
unknown measurement errors just mentioned. Note that these measurement 
errors are independent from scan to scan. That is, this variance will 
cause addition of white noise to the true time series* 

Since the variance of the true time series is duo to vary low 
frequency components (of long spatial scale), the true time series will 
fluctuate slowly and smoothly. This can be easily seen by noting that the 
10-km velocity integrations from one scan to. the next aro separated by 
only about 400-m horizontally. This separation is quite small when 
compared with the several-km scales contributing to the variance of the 
true time series. Clearly the added variance due to the independent 
errors will contribute short-term fluctuations to the time series, while 
the true time series will contribute predominantly long-term 
fluctuations. 


The solution to eliminating the error-contribution to the time 
series is obviously filtering. A low-pass filter wilL remove most of the 
contribution of the white errors, leaving most of the contribution of the 
true means. Note that this argument relies to a certain extent upon 
isotropy of the horizontal flow field, but only on scales smaller than l 
km. The most damage such an assumption could do to derived wind fields 
would be a reduction in the magnitude of l-kra and smaller flow features, 
and then the reduction would be suffered in only one dimension. 

The errors of this approach can easily be evaluated with sample wind 
fields. It holds promise of reducing velocity errors to very low levels, 
a requirement for processing boundary- layer fields for convective 
studies. 

3. Terrain returns. A careful examination of terrain returns will 
allow an independent check upon platform-velocity correction, since 
terrain would be expected to have zero velocity. Initial checks have 
shown terrain velocities on the order of 0-0.3 m/s. Terrain returns also 
have the potential for calibrating the lidar scanner both in elevation 
and azimuth. 

4. Anomalous errors. Unexplained errors remain in the data sets. 
These errors probably are due to hardware problems such as 
local-oscillator drift. These errors may be correlated with aircraft 
accelerometer data. Whatever the cause investigation is warranted. 


Appendix A: Magnetic tape format 


The following is a description of the format of the raw data tapes* 
Comments on decoding the header parameters are included. The tape records 
contain 352 16-bit words, written in DEC format. 


Word 

Name 

LSB 

Comments 

01 

X 

20 m 

East distance 

02 

Y 

20 m 

North distance 

03 

LOSD 

0.1 deg 

Add 180 deg 

04 

STATUS 

Coded 

See Integraph 81-029 

05 

PULSE WIDTH 

Coded 

See Integraph 81-029 

06 

PS 

Coded 

Processor status 

07-11 

SPARES 

NA 

NA 

12-299 

DATA 

- 

See text 

300 

NFL 

1 

Flight number 

301 

IDAY 

1 

Day number - Julian 

302 

NRUN 

1 . 

Run number 

303 

TSEC 

0.1 sec 

Seconds 0-59.9 

304 

TMIN 

1 min 

Minutes past midnight Z 

305 

PALT 

16.48 ft 

Pressure altitude 

306 

PALT 

7 

• 

Not significant 

307 

RALT 

2 ft 

Radar altitude 

308 

DFP 

0.1 deg C 

Dew/frost point 

309 . 

SATM 

0.1 deg C 

Static air temperature 

310 

IR-SUR 

0.1 deg C 

IR surface temperature 

311 

TAT 

0.1 deg C 

Total air temperature 

312 

LAT 

0.1 min 

Latitude 

313 

LAT 

1 deg 

Latitude 

314 

LON 

0.1 min 

Longitude 

315 

LON 

1 deg 

Longitude 

316 

TH 

0.04395 deg 

True heading 

317 

TAS 

0.5144 m/s 

True air speed 

318 

GS 

0.0628 m/s 

Ground speed actually used 

319 

DA 

0.04395 deg 

Drift angle 
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320 

WS 

0.05144 m/s 

INS wind magnitude 

321 

WD 

0.1 deg 

INS wind direction 

322 

P 

0.04395 deg 

Pitch 

323 

R 

0.04395 deg 

Roll 

324 

THETA 1 

0,1 deg 

Inner wedge positon 

325 

THETA2 

0.1 deg 

Outer wedge position 

326 

TTP 

10 ms 

Wedge time to position 

327 

SMT 

? 

Motor temperature 

328 

LOSE 

0.1 deg 

Scanner elevation 

329 

NINT 

1 

Pulses per integration 

330 

NLAG 

1 

Number of lags 

331 

TK 

0.04395- deg 

Track angle 

332 

LO 

0.01 MHz 

LO frequency 

333 

LOOFF 

0.08 m/s 

LO offset 

334 

GSSOURCE 

1 

GS source 0-2 

335 

GSD 

0.0628 m/s 

GS via Doppler radar 

336 

GSDA 

0.0628 m/s 

GS via Doppler/Addas 

337 

GS1A 

0.0628 m/s 

GS via INS/Addas 

338-351 

SPARES 

NA 

NA 

352 

GKSM 

1 

Checksum 

Note that 

DA, P, and 

R are 0-359 deg. They should be converted to 

bipolar values 

by subtracting 8192 if greater 

than 4096. The sa»e is true 

of LOSD and LOSE, from wh 

ich 3600 should be subtracted if they are 

greater than 

1800. See 

the comments in the 

program FCONVT in Appendix B 

for further conversion information. 


The data 

contained 

in words 12-299 

is composed of 96 sets of 3 


measurements, corresponding to the 96 320-m range gates. The first word 
of each set is a logarithmic amplitude (LSB=0. 184 dB) . The* second is a 
bipolar velocity (LSB*0.08 m/s). The third is a coded width estimate 
(0-15). 

TSGC and TMXN are subject to complex errors in flight 13. TK is not 
independent but is derived from DA and TH. 
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Appendix B: Data reduction software 


The following are brief descriptions of several programs developed 
during the course of the investigation of the data properties. While they 
constitute a complete set of data reduction software, it is not expected 
that they will be used as such. Rather, they may be of some use in 

helping to expain the technique* of data conversion, correction and 

* 

evaluation. All programs are written in DEC RT-11 Fortran. The plotting 
routines require a Hewlett-Packard 7221A plotter. 

1) Program FXLED1. This program allows the user to examine raw data 
files on tape and to plot any item at any scale. Optimum coefficients for 
drift-angle correction can be determined for any data set length, based 
upon a 4-lag correction filter. It is also possible to plot certain 
differential quantities. 

' 2) Program FTAPED. This program converts raw data files on tape to 
568-word disk files. These disk files can also be examined, dumped or 
plotted. 

3) Program FCONVT. This program operates on files produced by FTAPED 
and produces identical disk files as output with drift-angle corrections 
applied, X- and V-eoordinates corrected and added to the file, and mean 
powers and velocities computed. Plotting is also possible. 

4) Program FEDIT. FEDIT operates on files produced by FCONVT and 
creates a file identical except for the addition of standard-deviation 
estimates for each velocity estimate. 

5) Program FSMOTH. This program performs quadratic smoothing on 
files produced by FEDIT. 

6) Program FGRID. FGRID produces gridded flow-field plots or files, 
operating on files created by FSMOTH. It can also produce gridded plots 
from files created by FEDIT, with nearest-neighbor algorithms used 
instead of interpolation. 
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7) Subroutines GRAPH 1, SINV and MFSD. GRAPH! is a plotting 
subroutine used by most of the above programs. It ia designed to drive a 
Hewlett- Packard 7221A plotter. SINV and MFSD are used for matrix 
inversion in program FILED1. 



c 


PROGRAM FXLED1 

REVISED 10-22-81 FOR NEW DRIFT DEFINITION 


ORIGINAL PAGE W 
OF POOR QUALITY 


INTEGER*2 DBLK(4) , ERROR(4) , ARY ( 362 ) , CODE, NWDS , EKnADR 
INTEGER* 2 IRANGE ( 3 ) , IWID ( 3 ) , VELD ( 3 ) , IDRIFT ( 4 ) 

REAL* 4 AMP( 3 ) , VEL( 3 ) , WD(4) , B1 ( 3 ) , A1 (600) ,ATA(6) ,B(200) 
REAL*4 XI ( 3 ) , DRIFT1 ( 9 ) 

REAL* 8 TIME, FLEFT 
BYTE A 

DATA DBLK/3RMT0, 0,0,0/ 

DATA ERROR/4*0/ 

DATA WD/0.,1.,0.,-1./ 

ERRADR=I ACDR( ERROR ) 

ICHAN=IGETC() 

IF ( ICHAN.LT. 0) STOP 'NO CHANNEL AVAILABLE' 

IF ( LOOKUP ( ICHAN, DBLK, 0, -1 ) . LT . 0 ) STOP 'BAD LOOKUP' 

NRECD=0 

NSUM=103 

ZERO=-20. 

FULLS=20 . 

FMULT=. 00267 
SHIFT=0. 

ITAP13=0 
CONI 3 . 04395 
NWDS=352 
FACTO- 2 . 66 
LAST=1 
LAGV=3 

N RANGE. 3 3 6 . . 

10 TYPE 905 , NRECD 
ACCEPT * , I 


IF 

(1.EQ.0) 

GO 

TO 

90 

l READ AND TYPE HEADER 

IF 

(I .EQ.l) 

GO 

TO 

100 

1 READ A RECORD 

IF 

(I.EQ.2) 

GO 

TO 

200 

l PRINT A RECORD 

IF 

(I.EQ.3) 

GO 

TO 

300 

1 FORWARD SPACE 

IF 

(I.EQ.4) 

GO 

TO 

400 

l REWIND TAPE 

IF 

(I.EQ.5) 

GO 

TO 

500 

J BACKSPACE 

IF 

(I.EQ.6) 

GO 

TO 

600 

1 PLOT ROUTINE 

IF 

(I.EQ.7) 

GO 

TO 

700 

1 PLOT BOX 

IF 

(I.EQ.8) 

GO 

TO 

800 

1 SET WEIGHTS AND SCALE 

IF 

(I.EQ.9) 

GO 1 

TO 

350 

I SPACE TO RECORD N 

IF 

(I.EQ.10) 

GO 

TO 

801 

l SET LAGV, NSUM 

IF 

(I.EQ.11) 

GO 

TO 

750 

l OPTIMUM COEFS 

IF 

(I.EQ.12) 

GO 

TO 

802 

! SET # RANGES 

IF 

(I.EQ.13) 

GO 

TO 

803 

1 SET ITEM TO PLOT 

IF 

(I.EQ.14) 

GO 

TO 

804 

l SET PLOT TIME 

IF 

(I.EQ.15) 

GO 

TO 

905 

1 SET FMULT, SHIFT 

IF 

1I.EQ.16) 

GO 

TO 

806 

i SET PLOT SCALE 

IF 

(I.EQ.17) 

GO 

TO 

807 

l SET I TAP 13 - TIME CODE 

IF 

(I.EQ.18) 

GO 

TO 

750 

l PLOT OPTIMUM COEFS 

IF 

(I.EQ.19) 

GO 

TO 

750 

ITYPE OPTIMUM COEFS 

IF 

(I.EQ.99) 

GO 

TO 

999 

1 EXIT FROM PROGRAM 

GO 

TO 10 
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90 ITEST=0 I READ AND TYPE HEADER 

GO TO 105 
100 ITEST=1 
105 ERROR(l)=0 ' 

ERROR(2)“0 

ICODE=ISPFNW( 248 , ICHAN , NWDS , ARY, ERRADR) 

NRECD=NRECD+1 

ARY( 304 )=ARY( 304) .AND. 2047 
ARY( 303 )=ARY(303 ) .AND. 1023 
ITl=ARY(304)/60 
IT2=ARY( 304) -60*IT1 

T3=ARY( 303 )/l0 . I DECODE TIME 

IX=20*ARY(1) 

IY=20*ARY(2) IX, Y POSITIONS 

N1=ARY(313) 

N2=ARY ( 3 1 2 ) / 10 I LATITUDE 

IWl ss -ARY( 315) 

IW2=-ARY(314)/l0 1LONGTITUDE 
I FLT=ARY (300) 

I RUN-ARY (302) l FLIGHT AND RUN 

IALT=*ARY( 305 ) *16 . 4 l ALTITUDE 

rLOS*ARY(3)/l0.+180. 1 LINE OF SIGHT ANGLE 

THDG=ARY (316)* CONI l TRUE HEADING 
DRIFT=ARY(319)*C0N1 1 DRIFT ANGLE 
IF ( DRIFT. GT. 130. ) DRIFT=DRIFT-360 
FLOSE=ARY(328)/10. 1LOS ELEVATION 
PITCH=ARY( 322 ) *CONl t PITCH ANGLE 
ROLL= ARY (323)* CON 1 l ROLL ANGLE 
IF (ROLL. GT. 180. ) ROLL=ROLL-360. ' 

TAS=ARY( 317)*. 5144 1 TRUE AIRSPEED (M/S) 

GS=ARY( 318) *.0628 1 GROUND SPEED (M/S) 

WIND=ARY( 320)* .05144 1 INS WIND VELOCITY 

DI RN=ARY ( 3 2 1 ) / 1 0 . I INS WIND DIRECTION 
IPROC=ARY( 5 ) 1 PROCESSOR STATUS 

ISTAT=ARY ( 4 ) 1 MAJOR STATUS WORD 

ALOF= . 0 1 * ARY (332) 1LO FREQUENCY 
OFF^.GS^ARYOSS) ILO CORRECTION, M/S 
GS1=. 0628* ARY( 335 ) IDOPPLER GS VIA A/C 
GS2=.0628*ARY(336) IDOPPLER VIA ADDAS 
GS3=.0628*ARY(337) l INS GSPEED 
DRIFT2=CONl* (ARY ( 331 ) -ARY ( 316 ) ) 1 TRACK -TRHDG 

110 TYPE 901 , IX, IY, N1 , N2 , IW1 , IW2 , IFLT , IRUN, IT1, IT2 ,T3 , I ALT 
TYPE 902 , FLOS , THDG , DRIFT , FLOSE , PITCH , ROLL 
TYPE 90 3 , TAS , GS , WIND , DI RN , IP ROC , ISTAT 
TYPE 921 , ALOF, OFF, GS1 , GS2, GS3 , DRIFT2 
TYPE * 

IF (ITEST.EQ.0) GO TO 410 
I3=NRANGE/3 

DO 150 1=1,13 l PRINT NRANGE RANGE GATES 

DO 160 J= 1 , 3 

IRANGE( J )=320* 1+320*13 *J-3710 
I4=3*I3*J+3*I 



ORIGINAL PAGE IS 
OF POOR QUALITY 



AMP(J)=.134*ARY( 14-27) ORIGINAL PAGE IS 

VEL(J)=»08*ARY( 14-26 ) OF POOR QUALITY 

XWID(J)!=»ARY( 14-25) 

160 CONTINUE 

150 TYPE 904,(IRANGE(K),AMP(K),VEL(K),IWID(K),K=1,3) 

TYPE * 

GO TO 410 

200 PRINT 901,IX,IY,N1,N2,IW1,IW2,IFLT,IRUN,IT1,IT2,T3,IALT 
PRINT 902, FLOT,THDG, DRIFT, FCLOSE, PITCH, ROLL 
PRINT 903, TAS,GS, WIND, DIRN, IPROC,ISTAT 
PRINT * 

DO 250 1=1,12 
DO 260 J=1 , 3 

IRANGE(J)=320*I+3840*J-3710 
AMP (J)=. 184* ARY(36*J+3* 1-27) 

VEL(J)=.08*ARY(36*J+3*I-26) 

IWID(j)=ARY(36*J+3*I-25) 

260 CONTINUE 

250 PRINT 904, ( IRANGE(K) , AMP(K) ,VEL(K) , IWID(K) »K=1, 3) 

PRINT * 

GO TO 10 

300 TYPE 906 

ACCEPT * , N 
IF (N.EQ.0) GO TO 10 
310 ERROR ( 1 )=0 
ERROR ( 2 )=0 

ICODE=ISPFNW( 254 , ICHAN, N, ARY, ERRADR) 

NRECD=NRECD+N 

GO TO 410 ' • 

350 TYPE 917 1 SPACE TO RECORD N 

ACCEPT * , N 
N=N-NRECD 

IF (N.GT.0) GO TO 310 
IF (NnEQ.0) GO TO 10 
N=-N 

GO TO 510 

400 TYPE *, 'REWIND TAPE' 

ERROR ( 1 ) =0 
ERROR ( 2 ) =0 

ICODE=ISPFNW( 251, ICHAN, NWDS, ARY, ERRADR) 

NRECD=0 

410 IF (ICODE.NE.0) TYP E * , ' ICODE= ' , ICODE , ERROR ( 1 ) , ERROR ( 2 ) 
GO TO 10 

500 TYPE 907 

ACCEPT * , N 
IF (N.EQ.0) GO TO 10 
510 ERROR ( 1 ) =0 
ERROR ( 2 )=0 

ICODE=ISPFNW(253, ICHAN, N, ARY, ERRADR) 

NRECD=NRECD-N 
IF (NRECD.LT.0) NRECD=0 
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ORIGINAL PAGE IS 

IP (IC0DE.EQ.1) NRECD-0 OF POOR QUALITY 

GO TO 410 V ' T 

600 TYPE 914 l EXECUTE PLOT 

ACCEPT * , MATCH 1 MATCH THIS STATUS WORD (0»ANY) 

IF ( MATCH, LT.0) GO TO 10 

XFACTR=2000./(XAXIS*60. ) 1PUNITS PER SECOND*/- YFACTR=1400 . / ( FULLS—', 

A»'p' 

TYPE 915,155 l START PLOTTER \ 

DO 650 U»l, IRECDS 

ERROR( 1 ) — 0 
ERROR(2)=0 

ICODE=ISPFNW( 248 , ICHAN , NWDS , ARY, ERRADR) 

NRECD=NRECD+1 

ISUM=0 l SUM SEVERAL RANGES 

DO 610 K— 16, NSUM, 3 
ISUM=ISUM+ARY(K) 

610 CONTINUE 

ARY( 353 )=I5UM 

ARY( 303 )=ARY( 303) .AND. 1023 

ARY(304)=ARY(304) .AND. 2047 l MASK UNUSED BITS 

IF (ARY (319) .GT.4096) ARY( 319 )=ARY( 319 )-8192 
IF (ARY(322) .GT.4096) ARY(322)=ARY(322)-8192 
IF (ARY(323) .GT.4096) ARY( 323 )=ARY( 323) -8192 
IF (ARY( 328 ). GT . 1800 ) ARY ( 328 ) =ARY ( 328 ) -3600 

ARY(354)=ARY(335)-ARY(336) IA/C DOPPLER-ADDAS DOPPLER 

ARY (355) =ARY (335) -ARY (337) lA/C DOPPLER-ADDAS INS GS 
ARY ( 356 )=ARY( 336 ) -ARY( 337 ) IA/C ADDAS DOPPLER-ADDAS INS 
ARY ( 358 ) =ARY( ITEMD) -IDELY1 1 AND ANY OTHER ITEM 

ARY (359 )=ARY( ITEMD) -IDELY2 
ARY( 360 ) =ARY( 359 ) +ARY( 358 ) 

IDELY2=IDELY1 
IDELY1=ARY( ITEMD) 

IDRIFT ( 4 ) =IDRIFT ( 3 ) 

IDRIFT(3)=IDRIFT(2) 

IDRIFT ( 2 )=IDRIFT(1) 

IDRIFT ( 1 ) =ARY (331) -ARY ( 316) 1NEW DRIFT DEF'N - (ARY(319) IS OLD 

* DDRIFT=0. I WEIGHT DRIFT 

DO 630 1=1,4 

DDRIFT=DDRIFT+WD ( I ) * IDRIFT ( I ) 

630 CONTINUE 

ARY (357) = FACT R* DDR I FT *ARY ( 3 1 8 ) / 1 7 6 . 2 16+VELD ( LAGV ) 

VELD ( 5 )=0 .0 
VELD ( 4 ) =VELD ( 3 ) 

VELD ( 3 ) =VELD ( 2 ) 
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V ELD ( 2 ) »VELD ( l ) 
VELD(1)»ARY(353) 


ORIGINAL 
OF POOR QUALITY 


IF (MATCH. EQ.0) GOTO 640 
IF ( (ARY(4) .AND. 24) .NE. MATCH) GOTO 650 
640 TYPE 916 

TIME=*60 • D0*ARY( 304) + . 1D0*ARY( 303 ) 

IF (ITAP13.EQ.0) GO TO 647 ICODE TO FIX TAPE 13 
ISEC=ARY(303) / 10 . 

IF (MOD(lSEC,10) .EQ.9) TIME=TIME-1 . D0 
IF (MOD(ISEC,10) .NE.0) GO TO 645 
IF ( LAST . NE . 0 ) TIME=TIME-1.D0 
IF (LAST. ME. 0) ISEC-ISEC-1 
LAST=*0 
GO TO 646 

645 LAST=1 

646 CONTINUE 

IF (ISEC.LT.10) TIME=TIME+60.D0 

647 CONTINUE 

IF (TIME. LT. 43200. ) TIME«TIME+86400 . 

X=*XFACTR* ( TIME-FLEFT ) 

Y» ( ARY( NITEM ) * FMULT+SHIFT-ZERO ) ^YFACTR 
IF (X.GT.2000. ) X=2000. 

IF (Y.GT.I400.) Y=1400 . 

IF (Y.1-T.0. ) Y=0 . 

IF (X.LT.0.) A*'p' 

CALL GRAPH1(X,Y,A) ... . 

A=' q' 

650 CONTINUE 

TYPE 919,155 GO TO 10 

700 TYPE *, 'OUTLINE BOX' 

TYPE 915,155 

CALL GRAPH1 ( 0 . , 0 . , ' p ' ) 

CALL GRAPH1 ( 2000 . , 0 . , ' q ' ) 

CALL GRAPH1 ( 2000 . , 1400 . , ' q \ 

CALL GRAPH1(0., 1400., 'q') 

CALL GRAPH1 ( 0 . , 0 . , ' q ' ) 

TYPE 919,155 
GO TO 10 

750 ICASE=I 
ITRIAL=0 

TYPE * , ' FIND COEFS. FOR CODE (8,24):' 

ACCEPT * , MATCH 
IF (MATCH. LT.0) GO TO 10 
TYPE *, 'NUMBER OF SAMPLES:' 

ACCEPT * , NLSQ 
IF (ICASE.EQ.il) GO TO 754 
TYPE * , ' NUMBER OF TRIALS ' 

ACCEPT * , NTRIAL 

IF ( ICASE. EQ.19 ) GOTO 754 
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TYPE * , * PLOT COEFFICIENT 1-4 : ' 

ACCEPT * , NCOEF 
A*' p' 

TYPE 915,155 I TURN ON PLOTTER 

754 N— 9 

755 ERROR ( 1 ) =0 
ERROR( 2 )=0 
ICODE=ISPFNW( 248 , ICHAN, NWDS , ARY, ERRADR) 

NRECD=NRECD+1 

IF (ICODE.NE.0) GO TO 777 

IF (ARY (319) .GT.4096) ARY( 319 )=ARY( 319) -8192 
GS=.06279*ARY(318) IM/S UNITS 

DRIFT1 ( 9 ) =DRIFT1 ( 8 ) 

DRIFT1 (8 )=DRIFT1( 7 ) 

DRIFT! ( 7 )=DRIFT1 (6 ) 

DRIFT1(6)=DRIFT1(5) 

DRI FT 1 ( 5 ) =DRIFT 1(4) 

DRIFT1(4)=DRIFT1( 3) 

DRIFT1 ( 3 ) =DRIFT1 (2 ) 

DRIFT1 ( 2 )=DRIFT1( 1 ) 

DRIFT 1 ( 1 )= .04395* (ARY( 331 ) -ARY(316 ) ) *GS* .01645 IM/S UNITS 
C ABOVE LINE REVISED FROM (ARY (3 19)) IN PARENS 

IVEL9=IVEL8 

IVEL8=IVEL7 

IVEL7=IVEL6 

IVEL6=IVEL5 

■ IVEL5=IVEL4 • • ; •••*' • ••• • v • *- '*• • ' '' • • 

IVEL4=IVEL3 
IVEL3=IVEL2 
IVEL2=IVEL1 
IV EL 1=0 

DO 760 K=16,NSUM,3 
IVEL1 wsIVELl+ARY ( K) 

760 CONTINUE 1.00267 M/S UNITS 

IF (N.LT.l) GO TO 775 

IF ( MATCH. EQ.0) GO TO 765 
IF ( (ARY(4) .AND. 24) .NE. MATCH) GO. TO 755 
765 B(N)*'.00267*( 2*IVEL6-IVEL4-IVEL8 ) IM/S UNITS 

DO 770 J=1 , 3 
N1=N+NLSQ*(J-1) 

A1 ( N1 ) =DRIFT1 ( J ) -DRI FT 1 ( J+l ) -2 . *DRIFT1 ( J+2 ) +2 . *DRIFT1 ( J+3 ) 

A1 (N1 )=A1 (Nl)+DRIFTl(J+4) -DRIFT 1 ( J+5 ) 

770 CONTINUE 

775 N=N+1 

NP1=NLSQ+1 

IF (N.LT.NP1) GO TO 755 

CALL GTPRD( A1 , B, Bl, NLSQ, 3,1) 

CALL MATA(A1,ATA, NLSQ, 3,0) 
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ORIGINAL PAGE IS 

of poor quality 


CALL SINV(ATA,3, .0001, IER) 

CALL MPRD(ATA,B1,X1,3,3,1,0,1) 


ORIGINAL PAG* tt 
OF POOR QUALITY 


WD(1)=X1(1) 

VfD ( 2 ) “XI ( 2 ) -XI ( 1 ) 

WD ( 3 ) =*X1 (3)-Xl(2) 

WD(4)= S -X1 ( 3 ) 

ITRIAL»ITRIAL+1 

777 IF (ICODE.NE.0) ITRIAL=NTRIAL 

IF (ICASE.EQ.I1) GO TO 790 ICRT TYPEOUT 

IF (ICASE.EQ.19) GO TO 780 l PRINTOUT 

X»10.*ITRIAL 

Y® 700 . +WD ( NCOEF ) * 3 50 . 

IF (Y.GT.1400.) Y*1400. 

IF (Y.LT.0.) Y®0. 

CALL GRAPH1 (X, Y, A) 

A®' q' 

IF ( ITRIAL .LT .NTRIAL) GO TO 754 
TYPE 919,155 (TURN OFF PLOTTER 

GO TO 10 

780 ITl=ARY(304)/60 

IT2=ARY( 304)=60*IT1 
T3»ARY(303)/10. 

PRINT 920 , ITRIAL, NRECD, IT1 , IT2 , T3 , WD 
IF (ITRIAL.LT. NTRIAL) GO TO 754 
i, JO TO 10 

790 *1 YPE *■, 'WEIGHTING VECTOR*' , XI 

X (PE * , ' COEFFICIENTS® ' , WD 

GC TO 10 

800 TYl'E * , ' WEIGHTING VECTOR (.3 1 0 -1) ' 
ACCEPT * , WD 

TYPE *, 'SCALE FACTOR (3)' 

ACCEPT * , FACTR 
GO TO 10 

801 TYPE *, 'VELOCITY LAG (3)’ 

ACCEPT *,LAGV TYPE *,'NSUM (103)' 

ACCEPT * , NSUM 

TYPE * , ' ITEM TO DELAY' 

ACCEPT * , XTEMD 
GO TO 10 

802 TYPE *,' NUMBER OF RANGES (36)' 

ACCEPT * , NRANGE 

GO TO 10 

803 TYPE 908 

ACCEPT * , NITEM l ITEM TO PLOT 

GO TO 10 
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304 


805 


806 


807 


901 


902 

903 

904 

905 

906 

907 

908 

909 
9X0 

911 

912 

913 

914 

915 

916 

917 
913 

919 

920 

921 


999 


TYPE 911 

ACCEPT * , IHQUR, MINUTE 1LEFT EDGE INFO 
FLEFT-3600 . D0*IHOUR+60 . D0*MINUTE 
TYPE 912 

ACCEPT *,XAXIS IFS IN MINUTES 

TYPE 913 

ACCEPT * , IRECDS 1# RECORDS TO PLOT 

GO TO 10 

TYPE 910 

ACCEPT *,FMULT, SHIFT l PLOT FACTOR 

GO TO 10 

TYPE 909 

ACCEPT *, ZERO, FULLS 1 PLOT SCALE 
GO TO 10 

TYPE *, 'TAPE 13? (0,1) ' 

ACCEPT * , ITAP13 
GO TO 10 

FORMAT (IX, *X®' ,16, ' Y=',I6,' N * , 12 , • : ' , 12 , ' 

113,’: ',12,* FLT ',12,* RUN ’,13,’ TIME *,12, 
1* : ' ,12, ' : ' ,F4.1, ' ALT ' , 16) 

FORMAT (IX, 'LOS®' ,F6.1, ' THDG®' , F5 . 1 , ' DRIFT* 
1 ' LOSEs * , F5 . 1 , ' PITCH® * , F6 . 2 , ' ROLL® ' , F6 . 2 ] 
FORMAT (IX, 'TAS®' ,F5.1, ' GSPD® ' , F5 . 1 , ' WIND® 1 
1' WDIR® ' , F5 • 1 , ' PROC® ' , 06 , ' STATUS®', 06) 
FORMAT ( IX, 3 ( 15, F7 • 1 , F6 • 1, 14, 5X) ) 

FORMAT (IX, 'AT RECORD' ,15, ' COMMAND: ',$) 
FORMAT (IX, 'SPACE FORWARD N RECORDS: ',$) 

FORMAT (IX, 'SPACE BACKWARD N RECORDS: ',$) 
FORMAT (IX, 'PLOT ITEM M: ',$) 

FORMAT (IX, 'ZERO AND FULL-SCALE FOR PLOT: ',$) 
FORMAT (IX, 'MULTIPLIER AND OFFSET: ',$) 

FORMAT (IX, 'LEFT EDGE HOUR, MINUTE: ',$) 

FORMAT (IX, 'FULL SCALE IN MINUTES: ',$) 

FORMAT (IX, 'NUMBER OF RECORDS: ' ,$) 

FORMAT (IX, 'PLOT CODE TO MATCH (0,8,24): * ,$) 
FORMAT ( IX, A1 ,'•(',$) 

FORMAT (1X,A1,$) 

FORMAT (IX, 'SPACE TO RECORD N: ',$) 

FORMAT (IX, 5A1 , ' } ' ) 

FORMAT (IX, *p) ' ,A1, ' . ) ' ) 

FORMAT (IX, 14, 17, 15, ' : ' ,12, ' : ' , F4 . 1 , 4F10 . 4 ) 
FORMAT (IX, 'LO: ' ,F6.2, ' OFFSET: ' ,F6.1, ' GS:', 

1' DRIFT: ' , F5 . 2 ) 

CALL CLOSEC ( I CHAN ) 

STOP 'EXIT 99' 

END 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


w* , 

>' , F5 .2 , 
,F5.1, 


3F6.1, 





PROGRAM FTAPED 
REVISED 10-27-81 

TURKS MAGTAPE PILES INTO DISK PILES 
CAM ALSO EXAMINE/ PLOT TAPE DATA 


ORIGINAL Pfj® ® 
Of POOR QUALITY 


INTEGER* 2 DDLK(4) , ERROR (4 ) , ARY ( 354) , CODE, NWDS, ERRADR 
INTEGER* 2 I RANGE ( 3 ) , IWID( 3) 

REAL* 4 AMP (3 ) , VEL( 3) 

REAL* 8 TIME, FLEFT 
BYTE A, FILSPC( 11 ) 

DATA DBLK/3RMT0, 0,0,0/ 

DATA ERROR/ 4*0/ 

ERRADR»IADDR( ERROR) 

ICHAN«IGETC( ) 

IF ( I CHAN . LT . 0 ) STOP 'NO CHANNEL AVAILABLE' 

IF ( LOOKUP ( ICHAN, DBLK, 0, -1 ) .LT.0) STOP 'BAD LOOKUP* 

NRECD=»0 

FULLS=20. 

ZERO 8 *- 20 , 

LAST=*1 
CONI*. 04395 
FMULT*C0N1 
SHIFT®0. 

NWDS*352 

NRANGE=36 



ii 10 TYPE 905 , NRECD 
ACCEPT * , I 
IF (I .EQ.0) GO TO 90 
IF (I.EQ.1) GO TO 100 
IF (I.EQ.2) GO TO 200 
IF (I.EQ.3) GO TO 300 
IF (I.EQ.4) GO TO 400 
IF (I.EQ.5) GO TO 500 
IF (I.EQ.6) GO TO 600 
I IF (I.EQ.7) GO TO 700 

IF (I.EQ.8) GO TO 450 
IF (I.EQ.9) GO TO 350 
IF (I.EQ.12) GO TO 802 

IF (I.EQ.13) GO TO 803 

IF (I.EQ.14) GO TO 304 

IF (I.EQ.15) GO TO 805 

IF (I.EQ.16) GO TO 806 

IF (I.EQ.17) GO TO 807 

IF (I.EQ.99) GO TO 999 

GO TO 10 


1 READ AND TYPE HEADER 
l READ A RECORD 
1 PRINT A RECORD 
l FORWARD SPACE 
l REWIND TAPE 
l BACKSPACE 
1 PLOT ROUTINE 
1 PLOT BOX 
l CREATE A FILE 
l SPACE TO RECORD N 
l SET # RANGES 
1 SET ITEM TO PLOT 
l SET PLOT TIME 
l SET FMULT, SHIFT 
l SET PLOT SCALE 
l SET TAPE 113 - TIME CODE 
l EXIT FROM PROGRAM 


90 ITEST—0 l READ AND TYPE HEADER 

GO TO 105 
100 ITEST=1 
105 ERROR ( 1 )=0 
ERROR ( 2)= s 0 

ICODE=ISPFNW( 243 , ICHAN , NWDS , ARY, ERRADR) 
NRECD=NRECD+1 




4j 


IT1=ARY ( 304 ) / 60 
IT2=ARY(304)-60*IT1 
T3*ARY(303)/l0. 
IX=20*ARY( l) 
IY=20*ARY(2) 
N1=ARY(313) 
N2=ARY(312)/l0 
IW1=-ARY(315) 
IW2=-ARY( 314) /10 
IFLT=ARY(300) 
IRUN=ARY( 302 ) 
IALT=ARY( 305) *1.64 
FLOS=ARY(3)/l0.+180. 
THDG=ARY ( 316 ) *C0N1 
DRIFT=ARY( 319 ) *C0N1 
IF (DRIFT. GT. 130. ) 
FLOS E= ARY (3 2 8 ) / 1 0 . 


I DECODE TIME 
IX, Y POSITIONS 
I LATITUDE 
l LONGTITUDE 


l FLIGHT AND RUN 
1 ALTITUDE 

t LINE OF ■SIGHT ANGLE 
1TRUE HEADING 
1 DRIFT ANGLE 
DRIFT=DRIFT-360 
1 LOS ELEVATION 


PITCH=ARY ( 322 ) *CONl 1 PITCH ANGLE 
ROLL* ARY (323)* CON 1 l ROLL ANGLE 
IF ( ROLL. GT. 180. ) ROLL=ROLL-*360 . 
TAS=ARY(317)* .5144 1TRUE AIRSPEED (M/S) 
GS=ARY( 318)* .0628 1 GROUND SPEED (M/S) 

WIND=ARY( 320) *.05144 !INS WIND VELOCITY 


DI RN=ARY ( 3 2 1 ) / 1 0 . 

IPROC=ARY(5) 

ISTAT=ARY(4) 

ALOF= . 01*ARY( 332 ) 
OFF=.08*ARY(333) 
GS1=.0628*ARY(335) 
GS2= . 06 2 8* ARY ( 336 ) 
GS3 = .0628*ARY( 337) 


1 INS WIND DIRECTION 
! PROCESSOR STATUS 
l MAJOR STATUS WORD 
1 LO FREQUENCY 
ILO CORRECTION* M/S 
1 DOPPLER GS VIA- A/C- 
l DOPPLER VIA ADDAS 
l INS GSPEED 


DRIFT2=CONl* ( ARY ( 331 ) -ARY (316) ) ITRACK-TRHDG 


ORIGINAL PAGE IS 
OF POOR QUALITY 




110 TYPE 901* IX, IY,N1, N2, IW1* IW2 , IFLT, IRUN, IT1 , IT2 , T3 , IALT 
TYPE 902, FLOS, THDG, DRIFT, FLOSE , PITCH, ROLL 
TYPE 903 , TAS, GS , WIND, DIRN, IPROC, ISTAT 
TYPE 921, ALOF, OFF, GS1 , GS2 , GS3 , DRIFT2 
TYPE * 

IF (ITEST.EQ.O) GO TO 410- 
I3=NRANGE/3 

DO 150 1=1,13 1 PRINT NRANGE RANGE GATES 

DO 160 J=1 , 3 

I RANG E ( J )= 3 2 0 * 1+ 3 2 0 * 1 3 * J - 3 7 1 0 . ■ 

I4=3*I3*J+3*I 
AMP(J)=. 184* ARY (14-27) 

VEL( J)= .08*ARY( 14-26 ) 

IWID( J ) =ARY( 14-25 ) 

160 CONTINUE 

150 TYPE 904, (IRANGE(K) ,AMP(K) , VEL(K) , IWID(K) ,K=1, 3) 

TYPE * 

GO TO 410 

200 PRINT 901, IX, IY, N1 ,N2, IW1 , IW2, IFLT, IRUN, ITl , IT2 ,T3 , IALT 
* PRINT 902, FLOT, THDG, DRIFT, FCLOSE, PITCH, ROLL 
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PRINT 903, TAS,GS, WIND, DIRN, IPROC, ISTAT ■ _ bjw*c ia 

PRINT 92 1 , ALOF, OFF, GS1* GS2, GS3 , DRIFT2 ORIGINAL PAW K» 

PRINT * OF POOR QUALITY 

DO 250 1=1,12 
DO 260 J=1 , 3 

IRANGE(J)=320*It-3840*J-3710 
AMP ( J ) = . 184*ARY( 36*J+3*I-27 ) 

VEL(J)=.08*ARY( 36*0+3*1-26) 

IWID(j)=ARY(36*J+3*I-25) 

260 CONTINUE 

250 PRINT 904, ( IRANGE(K) , AMP(K) , VEL(K) , IWID(K) , K=l, 3 ) 

PRINT * 

GO TO 10 

300 TYPE 906 

ACCEPT *,N 
IF (N.EQ.0) GO TO 10 
310 ERROR ( 1 ) =0 
ERROR ( 2 ) =0 

ICODE=ISPFNW( 254 , ICHAN , N, ARY, ERRADR) 

NRECD=NRECD+N 
GO TO 410 

350 TYPE 917 1 SPACE TO RECORD N 

ACCEPT * , N 
N=N-NRECD 

IF (N.GT.0) GO TO 310 
IF (N.EQ.0) GO TO 10 
N=-N 

GO TO 510 .... - . .. V . . : • ... .• . .. ■ . 

400 TYPE “ f "REWIND TAPE' 

ERROR U}=0 
ERROR ( ^ ) =0 

lCQm~I S PFNW (251, ICHAN , NWDS , ARY , ERRADR ) 

NrECD=0 

410 IF (ICODE.NE.0) TYPE * , ' ICODE= ' , ICODE, ERROR ( 1 ji , ERROR ( 2 ) 

GO TO 10 

450 TYPE 451 I CREATE A NEW FILE 

451 FORMAT (IX, 'CREATE FILE DL0 : ' , $ ) 

ACCEPT 452, (FILSPC(I) , 1=1, 10) 

452 FORMAT (10A1) 

IF (FILSPC(l) .EQ.0) GO TO 450 
FILSPC ( 11 ) =0 
TYPE 453 

453 FORMAT (IX, 'NUMBER OF SCANS: ',?) 

ACCEPT * , NRECDS 

IF ( NRECDS. LE.0) GO TO 10 
N BLOCK*! . 3 94* FLOAT ( NRECDS ) + 1 
TYPE *, ' NBLOCK: ' , NBLOCK 

OPEN (UNIT=3, NAME=FILSPC, RECORDSIZE=177,TYPE='NEW' , 

.1 ERR=460 , FORM= ' UNFORMATTED ' , INITIALSIZE=NBLOCK) 

DO 455 1=1, NRECDS 
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error ( 1 ) “0 ORIGINAL PAGE IS 

ERROR( 2 ) =0 OF POOR QUALITY 

ICODE=ISPFNW( 243, ICIIAN, 352 , ARY, ERRADR) 

NRECD=NRECD+L j 

IF ( ICODE . ME . 0 ) GO TO 480 | 

ARY(353)=NRECD-1 f 

ARY( 354 )=0 !' 

WRITE (8) ARY 

455 CONTINUE i j 

V . 

456 CLOSE (UNIT=8) 

GO TO 10 


460 TYPE *, 'ERROR UPON OPENING' 

GO TO 10 

430 TYPE * , ' READ NON-STANDARD TAPE RECORD' 

GO TO 456 

500 TYPE 907 

ACCEPT * , N 
IF (N.EQ.0) GO TO 10 
510 ERROR ( 1 ) =0 ’ 

ERROR ( 2 ) =0 

ICODE=ISPFNW( 253, ICHAN,N, ARY, ERRADR) 

NRECD=N RECD-N 

IF ( Nl,„ CD . LT , 0 ) • ’NRECD=0 • • • ‘ *- ' ' ' / - V < 

IF (ISwDE.EQ.l) NRECD=0 
GO TO 410 

* ■ 

600 TYPE 914 l EXECUTE PLOT 

ACCEPT *, MATCH l MATCH THIS STATUS WORD (0=ANY) 

IF ( MATCH. LT.0) GO TO 10 

XFACTR=2000 . / (XAXIS*60 . ) ! PUNITS PER SECOND YFACTR=1400 . / ( FULLS-ZERO) 

A='p' 

TYPE 915,155 l START PLOTTER ’ 

DO 650 11=1, IRECDS 

ERROR ( 1 ) =0 
ERROR ( 2 ) =0 

ICODE=ISPFNW( 248 , ICHAN, NWDS , ARY, ERRADR) 

NRECD=NRECD+l 

IF (ARY (319) . GT. 4096) ARY ( 319) =- ARY ( 3 1 9 ) - 8 1 9 2 
IF (ARY (322). GT . 4096 ) ARY( 322 ) =ARY( 322) -3192 

IF ( ARY (323) .GT.4096) ARY( 323 ) =ARY( 323 ) -3192 ; ■>: '/ 

IF (ARY ( 323 ) .GT.1330) ARY ( 328 ) =ARY( 328 ) -3600 

IF ( MATCH. EQ.0) GOTO 640 • ( f, > 

IF ( ( ARY ( 4 ) . AND . 24 ) . NE . MATCH ) GOTO 650 * 

*640 TYPE 916 
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ORIGINAL PAGE is 

; : TIME»60 . D^^ARV‘(30A) + .1D0*ARY(3O3) OF POOR Q U * L 

IP ( ITAPl’3 . EQ . 0 ) GO TO 64j'7 l CODE TO FIX TAPE 13 
ISSC=ARY(303)/10, 1/ 

IF ( MOD ( I SEC, 10) .EQ.9) TIME*TIME-1 . D0 

IF (MOD(ISEC,10) .NE.0) GO TO 645 

IF (LAST. NE.0) TIME=TIME-1 . D0 

IF (LAST. ME. 0) ISECaISEC-1 

LAST®0 

GO TO 646 

645 LAST® I 

646 CONTINUE 

IF (ISEC.LT.10) TIME=TIME+60.D0 

647 CONTINUE 

IF (TIME. LT. 43200. ) TIME=TIME+86400 . y 

X=XFACTR* (TIME-FLEFT) 

Y= ( ARY ( N ITEM ) * FMULT+SHI FT— ZERO ) * YFACTR 
IF (X.GT.2000.) X=2000. 

IF (Y.GT.1400.) Y=1400. 

IF (Y.LT.0.) Y=0. 
j IF (X.LT.0.) A=’p' 

i 

f CALL GRAPE11 (X, Y» A ) 

\ I q * 

! 650 CONTINUE 

j TYPE 919,155 GO TO 10 

! 700 TYPE * , ' OUTLINE BOX' 

TYPE 915,155 ’ *• • - • * • •• * • • 

CALL 'GRAPH! ( 0. ,0. , 'p' ) 

CALL GRAPH! (2000 . , 0 . , ' q ' ) 

CALL GRAPH1 (2000., 1400., 'q* ) 

CALL GRAPH! (0., 1403., 'q* ) 

CALL GRAPH1 ( 0 . , 0 . , ' q ' ) - 

TYPE 919,155 

GO TO 10 *' 

802 TYPE NUMBER OF RANGES (36)' 

ACCEPT * , NRANGE 

GO TO 10 

803 TYPE 908 

ACCEPT * , NIT EM l ITEM TO PLOT 

GO TO 10 

904 TYPE 911 

ACCEPT *, IHOUR.MINUTF, l LEFT EDGE INFO 
FLEFT=3600 . D0*IHOUR+4>0 . D0* MINUTE 
TYPE 912 

ACCEPT * , XAXIS IFS IN MINUTES 

TYPE 913 

Vi ACCEPT *, IRECOS 1 f RECORDS TO PLOT 

GO TO 10 
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8.35 type 913 '{ Origin a i " 

ACCEPT * , FMULT , SHIFT 1 PLOT FACTOR ' OF P AG£ Iff 

GO TO L0 P0 °ft QUALITY 

806 TYPE 909 

ACCEPT *,2ERO, FULLS l PLOT SCALE 
GO TO 10 


807 TYPE * , ' TAPE 13? (0,1)' 
ACCEPT *, IT API 3 
GO TO 10 


931 


902 

903 

904 

905 

906 

907 

908 
939 

910 

911 

912 

913 
'914 

915 

916 

917 

918 

919 

920 

921 


FORMAT (IX, *X=' ,16, ' Y=',I6,' N ' , 12 , ' : ' , 12 , ' W , 

113,’:', 12,' FLT ',12,' RUN M3,' TIME ',12, 

1’ ' ,12, ' : ' ,F4.1, ' ALT ',I5,'0') 

FORMAT (IX, ' LOS* ' , F6.1, * THDG* ' , F5 . 1 , ' DRIFT* ' , F5 . 2 , 
1' LOSE*' , F5 .1 , ' PITCH* ' , F6 . 2 , ' ROLL*',F6.2) 

FORMAT (IX, 'TAS*' , F5.1, ' GSPD*' ,F5.1, ' WIND* ' , F5 . 1 , 
1’ WDIR*' , F5 . 1 , ' PROC* ' , 06 , ' STATUS*', 06) 

FORMAT (IX, 3(15, F7 . 1, F6 . 1, 14, 5X) ) 

FORMAT ( IX, ' AT RECORD ' , 15 , ' COMMAND : ' , $ ) 

FORMAT (IX, 'SPACE FORWARD N RECORDS: ',$) 

FORMAT (IX, 'SPACE BACKWARD N RECORDS: ',$) 

FORMAT (IX, 'PLOT ITEM M: ',$) 

FORMAT (IX, 'ZERO AND FULL-SCALE FOR PLOT: ',$) 

FORMAT (IX, 'MULTIPLIER AND OFFSET: ',$) 

FORMAT (IX, 'LEFT EDGE HOUR, MINUTE: ',$) 

FORMAT (IX, 'FULL SCALE IN MINUTES: ',$) 

FORMAT (IX, 'NUMBER OF RECORDS: ',?) 

FORMAT (IX, 'PLOT CODE TO MATCH (0,8,24): ',$) 

FORMAT (IX, Al, ' . ( ' ,?) 

FORMAT ( IX, Al , $ ) 

FORMAT (IX, 'SPACE TO RECORD N: *,$) 

FORMAT (IX, 5A1, ' } ' ) 

FORMAT (IX, * p) ' ,A1, ' .) ' ) 

FORMAT (IX, 14, 17, 15, ' : ' ,12, ' : ' » F4 . 1 , 4F10.4 ) 

FORMAT (IX, 'LO: ' ,F6.2, ' OFFSET: ', F6 . 1 , ' GS:',3F6.1, 

1* DRIFT: ' , F5 . 2 ) 


999 CALL CLOSEC ( ICHAN ) 
STOP 'EXIT 99' 

END 


Q O O O O O O 


PROGRAM. FCONVT . 


ORIGINAL PAGE IS 
OF POOR QUALITY 


(i 

if. 


REVISED 10-29-81, 11-5,16-01, 4-13-82 
OPERATES ON FTAPED FILES (177 WORD RECORDS) 
TO PRODUCE .RAW FILE^ (560 WORD RECORDS) 


CAN ALSO EXAMINE 568-WORD RECORDS AND PLOT THEM 


INTEGER* 2 HEADER( 56 ) , DATA( 12, 90 ) 

INTEGER*2 RAW (354) , RAW1 (354) , RAW2 ( 354) , RAW3 ( 354 ) 

REAL*8 TIME, FLEFT 

REAL*4 COEF1 , COEF2 , COEF3 , COEF4 

BYTE A, F1LSPC( 11) 

ITAP13=0 

FULLS=20. 

ZERO=-20. 

LAST-1 
FMULT=.01 
SHIFT-0 . 

, MRANGE=1 
NRANGE=36 

COEF1=0.32 1 DRIFT CORRECTION COEFFICIENTS 

COEF2-0 . 56 
COEF3=0 • 12 
COEF4=- 1 . 0 

10 TYPE 905 , NRECD 
ACCEPT * , I 



IF 

(I.EQ.0) 

GO TO 

90 

l READ AND TYPE HEADER 

■ 1 • . 

IF 

(I , EQ . 1 ) 

GO TO 

100 

l READ A . RECORD' \ . 


IF 

(I.EQ.2) 

GO TO 

200 

1 PRINT RECORD JUST READ 


IF 

( I . EQ . 4 ) 

GO TO 

300 

1 REWIND FILE 


IF 

(I.EQ.6) 

GO TO 

600 

1 PLOT ROUTINE 


IF 

(I.EQ.7) 

GO TO 

700 

l PLOT BOX 


IF 

( I . EQ . 8 ) 

GO TO 

400 

l CREATE A NEW FILE 


IF 

( I . EQ . 9 ) 

GO TO 

350 

1 SPACE TO RECORD N 


IF 

(I.EQ.10) 

GO 

TO 

750 

10PEN READ FILE 


IF 

(I.EQ.ll) 

GO 

TO 

770 

1 CLOSE READ FILE 


IF 

(I.EQ.12) 

GO 

TO 

802 

l SET # RANGES 


IF 

(I.EQ.13) 

GO 

TO 

303 

1 SET ITEM TO PLOT 


IF 

(I.EQ.14) 

GO 

TO 

804 

1 SET PLOT TIME 


IF 

(I .EQ.15) 

GO 

TO 

805 

ISET FMULT, SHIFT 


IF 

(I .EQ. 16) 

GO 

TO 

806 

1 SET PLOT SCALE 


IF 

(I.EQ.17) 

GO 

TO 

307 

ISET TAPE 113 - TIME CODE 


IF 

(I.EQ.18) 

GO 

TO 

808 

ISET COEFFICIENTS 


IF 

(I.EQ.99) 

GO 

TO 

999 

i EXIT FROM PROGRAM 


GO 

TO 10 





90 

ITEST=0 



1 READ 

i AND TYPE HEADER 


GO 

TO 105 





100 

ITEST=1 





105 

READ ( 0 , ERR= 

670 

, END=660) 

HEADER, DATA 


NRECD=NRECD+1 

ITl=HEADER(4)/60 
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IT2*HEADER(4)-60*IT1 

T3*HEALER(5)/l3* 

XX* 2* HEADER (6 ) 

IY-2*HEADER(7) 

N1*HEADER(8) 

N2*HEADER(9)/l0 

IW1=HEADER(10) 

IW2=HEADER(ll)/l0 

IFLT=HEADER(2) 

IRUN=HEADER(3) 

IALT=HEADER(12 ) *2 

FLOS*HEADER(31)/l0. 


1 DECODE TIME 
IX,Y POSITIONS 
1 LATITUDE 


ORIGINAL PAGE 19 
OE POOR QUALITY 


l LONGTITUDE 


RUN 


l FLIGHT AND 
l ALTITUDE 

1 LINE OP SIGHT ANGLE 


THDG=HEADER(26 ) * ,04395 
DRIFT=HEADER( 27 ) * . 04395 
FLOSE=IIEADER ( 3 2 ) / 1 0 . 

PITCH=HEADER( 29 ) * . 04395 
ROLL=HEADER( 30) *.04395 
TAS=HEADER( 20 ) / 100 . I TRUE 
G S*HE ADE R ( 42 ) / 1 0 0 . 

WIND=HEADER( 43 ) /100 . 

DIRN=HEADER(44)/10. 

IPROC=HEADER( 19 ) 

ISTAT=HEADER( 18 ) 

ALOF=HEADER( 3 3 ) / 1 00 , 

OFP=HEADER( 3 4 ) / 1 0 0 . 

GS1=HEADER( 21 ) / 100 . 

GS2=HEADER(22 ) /100 . 

GS3 s =HEADER(23 ) /100 . 

DRIFT2=* . 04395 * ( HEADER ( 28 ) -HEADER ( 26 ) ) 
FLOSA=HEADER ( 46 ) / 10 . 1 ACTUAL L.O.S, 

POWER=HEADER(36)/100. 1 MEAN NOISE POWER 

VAVG=HEADER(41 ) /100 . IMEAN VELOCITY 

VCOR=HEADER(40)/100. 1 VELOCITY CORRECTION 


l TRUE HEADING 
1 DRIFT ANGLE 
l LOS ELEVATION 
1 PITCH ANGLE 
l ROLL ANGLE 
AIRSPEED (M/S) 
l GROUND SPEED (M/S) 

l INS WIND VELOCITY 
HNS WIND DIRECTION 
l PROCESSOR STATUS 
l MAJOR STATUS WORD 
ILO FREQUENCY 
1 LO CORRECTION, M/S 
l DOPPLER GS VIA A/C 
l DOPPLER VIA ADDAS 
1 INS GSPEED 


1 TRACK-TRHDG 


110 TYPE 901,IX,IY,N1,N2,IW1,IW2,IFLT,IRUN,IT1,IT2,T3,IALT 
TYPE 902 , FLOS , THDG, DRIFT, FLOSE, PITCH, ROLL 

TYPE 903, TAS,GS, WIND, DIRN, IPROC, ISTAT 
TYPE 921, ALOF , OFF , GS1 , GS2 , GS3 , DRIFT2 
TYPE 922, FLOSA, POWER, VAVG , VCOR, HEADER( 1 ) 

TYPE * 

IF (ITEST.NE.2) GO TO 111 

PRINT 901 , IX, IY,N1 , N2 , IW1 , IW2, IFLT , IRUN , IT1 , IT2 ,T3 , IALT 

PRINT 902, FLOS , THDG , DRI FT , FLOSE , PITCH, ROLL 

PRINT 903, TAS,GS, WIND, DIRN, IPROC, ISTAT 

PRINT 921, ALOF , OFF , GS 1 , GS2 , GS3 , DRIFT2 

PRINT 922 , FLOSA, POWER, VAVG , VCOR, HEADER ( 1 ) 

PRINT * 

111 IF (ITEST.EQ.0) GO TO 10 

DO 150 J=MRANGE , NRANGE ! PRINT NRANGE RANGE GATES 

IRANGE=320*J+32. 

AMP=.01*DATA(1, J) 

VEL= .01* DATA ( 2 , J ) 

IWID=DATA(3,J) 
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, --a--: -fs— -T— : ' “rr ■ r. *4 .. .’•■ ■ ***** 


ORIGIN At PAGE IJ 
OF POOR QUALITY 


IX*2.*DATA(4,»J) 

• ’ IY*2. *DATA(5, J) 

**' ^^-^*DATA(7 > J) 

V2=.01*DATA(8,J) 

V3=.01*DATA(9, J) 

V4«.01*DATA( 10, J) 

V5= . 01 * DATA ( 1 1 , J ) 

SIG*,01*DATA(6,J) 

SIG2=,01*DATA(12,J) 

IF (ITEST.NE.2) GO TO 150 

PRINT 904# IRANGE, AMP,VEL, IWID, IX# IY, SIG, Vl , V2 , V3 , V4 , V5 , SIG2 
150 TYPE 904, IRANGE, AMp, VEL, IWID, IX, IY, SIG , VI, V2 , V3 , V4, V5 , SIG2 
TYPE * 

IF (ITEST.EQ.2) PRINT ♦ 

IF (ITEST.EQ.2) REVfIND 6 
GO TO 10 


l EDITING STD DEVN 
t SMOOTH STD DEVN 


•** 


* 


200 ITEST=2 l PRINT RECORD 

GO TO 105 

350 TYPE 917 l SPACE TO RECORD N 
ACCEPT * , N 

IF (N.LT.l) GO TO 350 
IF (N.GT.MAXRCD) GO TO 350 
IF (N.EQ.NRECD) GO TO 10 
IF (N.GT.NRECD) GO TO 360 
N1=NRECD-N 
DO 355 J=1,N1 
355 BACKSPACE 3 

.. . .GO. TO 370- • • ' 

360 N1=N-NRECD 

DO 365 J=1,N1 

365 READ ( 8 , ERR=670 , END=660 ) 

370 NRECD-N ; * 

GO TO 10 

300 TYPE *, 'REWIND FILE' 

REWIND 8 
NRECD=1 
GO TO 10 

400 TYPE 401 l CREATE A NEW FILE 

401 FORMAT (IX, 'CREATE FILE DL0:',$) 

ACCEPT 402, (FILSPC(I) , 1=1, 10) 

402 FORMAT (10A1) 

IF (FILSPC(l) .EQ.0) GO TO 400 
FILSPC (11) =0 
TYPE 403 

403 FORMAT ( IX, ' NUMBER* OF SCANS: ',$) 

ACCEPT * , NRECDS 

IF ( NRECDS. LE.0) GO TO 10 
NBLOCK=4 .473 * FLOAT ( NRECDS ) + 1 
TYPE * , ' NBLOCK: ' , NBLOCK 

OPEN ( UNIT=9 , NAME=FILSPC, RECORDSIZE=568 , TYPE= ' NEW' , 
1 ERR=460 , FORM® ' UNFORMATTED ' , INITIALSI ZE=NBLOCK ) 
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OF o™ L PAGE 13 
Of POOR QUAUJY 


TYPE 404 rwu " QUALM 

FORMAT (IX,* INPUT FILE; DL3 : ' , $ ) 

ACCEPT 402, (FILSPC(I), 1-1,10) 

FILSPC(ll)=0 

OPEN ( UNITES, NAME-FILSPC, RECORDSIZE-177 , TYPE- ' OLD' 
1ERR-470, FORM- ' UNFORMATTED * , READONLY) 

READ (8, ERR-480, END-465) RAW1 
READ ( 3, ERR-480, END-465) RAW 2 
READ (8, ERR-480, END-465) RAW 3 

DO 455 I- I, NRECDS 

DO 410 J-1,354 
RAW(J)-RAW1(J) 

RAW1 ( J ) -RAW2 ( J ) 

RAW2 ( J ) -RAW3 ( J ) 

READ (3, ERR-480, END-465) RAW 3 


CONVERSION GOES HERE 

HEADER ( l ) -RAW ( 353 ) 
HEADER (2 ) -RAW ( 300 ) 
HEADER ( 3 ) -RAW ( 302 ) 
HEADER ( 4 ) -RAW ( 304 ) 
HEADER ( 5 ) =RAW( 303 ) 
HEADER ( 6 ) -RAW ( 1 ) 
HEADER( 7 )-RAW( 2 ) 
HEADER(8 ) =RAW( 313 ) 
HEADER( 9 )=RAw! 312 ) 
HEADER (10) —RAW (315) 
HEADER(ll)— RAW (3 14) 
IF ( IABS ( RAW ( 305 ) ) .GT 
HEADER(12)=8. 2 35 7* RAW 
HEADER! 13 )=RAW( 307) 
HEADER! 14 )=RAW! 303) 
HEADER ! 1 5 ) -RAW ! 309 ) 
HEADER( 16 )— RAW( 310) 
HEADER(17 )— RAW(311) 
HEADER( 18 )— RAW(4 ) 
HEADER! 19) -RAW! 6) 
HEADER ( 20 )=51 . 44* RAW ( 
IF ( IABS ( RAW( 335 ) ) .GT 
HEADER(21)= 6.28*RAW(3 
IF ( IABS (RAW( 336 )) .GT 
HEADER ( 22 )=6 . 28*RAW( 3 
IF (IABS (RAW (337 ) ) .GT 
HEADER ( 23 )=6 . 28* RAW ( 3 
HEADER! 24 )=RAW( 324) 
IIEADER( 25 )=RAW( 325 ) 
HEADER(26 )— RAW( 316 ) 
HEADER ( 27 )=RAW!319) 
HEADER! 28) -RAW! 331 ) 


t RECORD # 
l FLIGHT * 
l RUN # 
l TIME, MIN 
l TIME, .1 SEC 
. 1X0, 20 M 
l Y0, 20 M 
1LAT, DEG 
lLAT, .1 MIN 

l LONG, DEG 
1 LONG, .1 MIN 
.3000.) RAW ( 305 ) =0 
(305) 1PALT, 2 FT 

l RALT , 2 FT 
l DEW FT, .1 DEG 
1TEMP, .1 DEG 
HR TEMP, .1 DEG 
l TOT AIR T, .1 DEG 
l STATUS WORD 
l PROCESSOR STATUS 
317) ITAS, .01 M/S 

.4500) RAW ( 3 3 5 ) -0 

35) IGS-DOPPLER, .01 M/S 
.4500) RAW( 336 )— 0 

36) IGS-DOP/ADDAS, .01 M/S 
.4500) RAW (337) =0 

37) l GS- INS/ ADD AS , .01 M/S 

l WEDGE (I), .1 DEG 

l WEDGE jo), .1 DEG 
iTHDG, .04395 DEG 
(DRIFT, .04395 DEG 
(TRACK, .04395 DEG 
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j|tfafaitf- -r- — <tgpzgr~ ;:~:-j liWf'- 


..x 


READER (29) 
HEADER( 30) 
HEADER ( 31 ) 
HEADER(32) 
HEADER(33) 
HEADER (34) 
HEADER( 35 ) 
HEADER (36) 
HEADER ( 37 ) 
HEADER(3Q) 
HEADER(39) 
HEADER(40) 
HEADER(41) 
HEADER(42) 
HEADER (43) 
HEADER(44) 
HEADER (45) 
HEADER (46) 
HEADER (47) 
HEADER (43) 
HEADER (49) 
HEADER (50) 
HEADER (51) 


“ RAW (322) 
«RAW(323) 

“RAW ( 3) +1800 
“RAW (328) 

“RAW ( 332 ) 
=8*RAW(333) 
“RAW( 306 ) 

=0 
“0 
=0 
=0 
= 0 
“0 

- 6 . 28*RAW( 318 ) 
=5 . 144* RAW ( 320 ) 
“RAW(321) 

“RAW ( 326 ) 

=0 

=0 

=0 

=0 

=0 

=0 


ORIGINAL PAGE IS 
OF POOR QUALITY 

l PITCH, ,04395 DEG 
l ROLL, .04395 DEG 

l LOS DIRN, .1 DEG 
l LOS ELEV, ,1 DEG 
l L • 0 . , .01 MHZ 

iL.O. OFST, .01 M/S 
IPALT, FINE 
l MEAN POWER, .01 DB 
l DRIFT CORN, .01 M/S 
1SMOTH CORN, .01 M/S 
l BAD SCAN FLAG 
IV CORN ADDED, .01 M/S 
i MEAN V, .01 M/S 
1GS USED, .01 M/S 

1INS WIND, .01 M/S 
i INS DIRN, .1 DEG 
l TIME TO POSN, .01 SEC 
1 ACTUAL LOS DIRN, .1 DEG 
l SPARES FOLLOW 


HEADER ( 52 ) =0 
HEADER ( 53 )“0 
HEADER ( 54) =0 ' 
HEADER ( 55 )=0 
HEADER (56 ) =0 


IF ( HEADER ( 27 ) .GT . 4096 ) 
IF (HEADER(29) .GT.4096) 
IF (HEADER( 30) .GT.4096) 
I F ( HEADER ( 3 1 ) . GT . 1 800 ) 
IF (HEADER(32) .GT.1800) 


HEADER( 27 )=HEADER( 27 ) -8192 
HEADER ( 29 ) =HEADER( 29 ) -8192 
HEADER ( 3 0 ) “HEADER ( 3 0 ) -8 1 92 
HEADER ( 3 1 ) “HEADER (31) -3600 
HEADER( 32 ) =HEADER( 32 ) -3600 


♦ 


POWER=0 . 

VEL=0 . 

DO 415 J=1 , 20 

POWER=POWER+FLOAT ( RAW (3*J+219))*18.4 
VEL“VEL+FLOAT ( RAW ( 3*J+1 3 ) ) *8 . 

415 CONTINUE 

POWER=POWER/20. 

IF (ABS(POWER) .GT. 20000. ) POWER=20000 . 

IF (ABS(VEL) .GT. 160000. ) VEL=0. 

HEADER (36) “POWER 
HEADER(41)=(VEL/20. ) 

ALOS=.l*FLOAT(RAW(3) )+180. t LOS DIRECTION, DEG 

THDGC=FLOAT ( RAW1 { 316 ) -RAW ( 316 ) ) ITHDG CHANGE, .04395 DEG 

IF ( THDGC .GT.4096.) THDGC“THDGC-8192 . 

IF ( THDGC . LT . -4096 . ) THDGC=THDGC+S1 9 2 . 
ALOS=ALOS+.75*.04395*THDGC l ACTUAL LOS DIRN, DEG 

CCOS=16 . *COS( ALOS/57 . 2958 ) IY INCREMENT PER RANGE, 20 M 

CSIN=16 . *SIN( ALOS/57 . 2958 ) IX INCREMENT PER RANGE, 20 M 

X0“ • l*CSIN+FLOAT ( RAW ( 1 ) ) IX FOR J=0 
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Y0=.1*CCOS+FLOAT(RAW(2) ) IY FOR J=0 ^OOf? QUALITY 

,, HEADER (46) =10. *ALOS 

Vl«COEFl*FLOAT ( RAW3 ( 319 ) ) l DRIFT CORRECTION 

V2»COEF2*FLOAT(RAW2(319) ) 

V3=COEF3 * FLOAT ( RAW1 ( 3 19 ) ) 

V4=COEF4* FLOAT ( RAW (319)) 

VCORN=FLOAT( HEADER (42) )/l383. 1 CORRECTION TO ADD TO V, .01 M/S 
VCORN=VCORN* ( V1+V2+V3+V4 ) 

IF ( ABS ( VCORN ) .GT. 2000. ) VCORN=0. I CORRECT IOVFL 

ICORN=VCORN 

HEADER( 37 )=ICORN 

HEADER ( 40 )=ICORN 

HEADER ( 41 ) =HEADER ( 41 ) +ICORN 

DO 420 J=l,12 
DO 420 K=l,90 
420 DATA ( J , K ) =0 

DO 430 J=l, 90 

TEMP=18 . 4* FLOAT ( RAW( 3*J+9 ) ) l AMP, .01 DB 
IF ( ABS (TEMP) .GT. 20000. ) TEMP=20000. 

DATA ( 1 , J ) =TEMP 

TEMP=8 . *RAW( 3*J+10 ) +ICORN 1VEL, .01 M/S 

IF ( ABS ( TEMP ) .GT. 8000. ) TEMP=0. 

DATA ( 2 ? J ) sTEMP 

DATA(3,J)=RAW(3*J+11) I WIDTH, CODED 

DATA(4, J)=X0+J*CSIN l EAST DISTANCE 

DATA ( 5 ,, J ) =Y0+ J * CCOS l NORTH DISTANCE 

430 CONTINUE.. . . , . . .. .. . . •. • . -V 

WRITE (9,ERR=490) HEADER, DATA 

455 CONTINUE 

456 CLOSE (UNIT=8,ERR=485) 

CLOSE (UNIT=9,ERR=485) 

GO TO 10 

460 TYPE *, ‘ERROR IN OPENING OUTPUT FILE' 

GO TO 10 

465 TYPE * , ' END OF FILE' 

GO TO 456 

470 TYPE *, 'ERROR IN OPENING INPUT FILE' 

GO TO 456 

480 TYPE * , ' READ ERROR' 

GO TO 456 

435 TYPE *, 'CLOSURE ERROR' • 

GO TO 10 

490 TYPE *, 'WRITE ERROR* 
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*>* 

T**- 

600 


GO TO 456 
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TYPE 914 I EXECUTE PLOT 

ACCEPT * , MATCH l MATCH THIS STATUS WORD ( 0= ANY ) 

IF' ( MATCH. LT.0) GO TO 10 


ft 


XFACTR=2000 . / (XAXI§*60 . ) l PUNITS PER SECOND YFACTR=1400 . / ( PULLS 
A= * p * 

TYPE 915,155 I ST ART PLOTTER 

DO 650 11=1, IRECDS 


READ (8, END=660,ERR=670) HEADER, DATA 
NRECD=NRECD+1 


640 


IP (MATCH. EQ.0) GOTO 640 

IP ( (HEADER(18) .AND. 24) .NE. MATCH) GOTO 650 
TYPE 916 

TIME=60 . D0*HEADER(4) + . 1D0*HEADER( 5 ) 


IF (ITAP13.EQ.0) GO TO 647 l CODE TO FIX TAPE 13 
ISEC= . 1*HEADER( 5 ) 

IF (M0D(ISEC, 10) .EQ.9) TIME=TIME-1 . D0 

IF ( MOD( ISEC, 10 ) . NE . 0 ) GO TO 645 

IF ( LAST . NE . 0 ) TIME=TIME-1.D0 

IF ( LAST . NE .0 ) ISEC=ISEC-1 

LAST=0 

GO TO 646 

645 LAST=1 

646 - CONTINUE . • - 

• IF (ISEC .LT .10) T T ME =T I ME +6 0*. D 0 • 

647 CONTINUE 

IF (TIME. LT. 43200. ) TIME=TIME+86400 . 

X=XFACTR* (TIME-FLEFT) 

Y= ( HEADER ( NITEM ) *FMULT+SHlFT-ZERO) *YFACTR 
IF (X.GT.2000.) X=2000. 

IF (Y.GT.1400.) Y=1400. 

IF (Y.LT.0. ) Y=0 . 

IF (X.LT.0. ) A='p* 

CALL GRAPH1(X,Y,A) 

A= 1 q 1 

650 CONTINUE 

TYPE 919,155 GO TO 10 

660 TYPE * , 'VEND OF FILE* 

GO TO 10 

673 TYPE * , ' READ ERROR' 

GO TO 10 

700 TYPE OUTLINE BOX' 

TYPE 915,155 

CALL GRAPH 1 ( 0 * , 0 . , ' p ' ) 

CALL GRAPH 1 ( 2000 . , 0 . , ' q ' ) 
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CALL GRAPHl (2000., 1400,, 'q' ) 0F P °OR QUALITY 

CALL GRAPHl (0., 1400., 'q' ) 

CALL GRAPH 1 (0 . , 0 . , ' q' ) 

TYPE 919,155 
GO TO 10 

750 TYPE 751 

751 FORMAT (IX, 'OPEN FOR READING FILE DL0 t ' , $ ) 

ACCEPT 752, (FILSPCd) , 1*1, 10) 

752 FORMAT (10A1) 

IF (FlLSPC(l) .EQ.0) GO TO 750 

FILSPC( 11 ) *0 

TYPE TOTAL RECORDS:' 

ACCEPT * , MAXRCD 

OPEN (UNIT*8 , NAME'fFILSPC, RECORDSIZE=568 , TYPE* * OLD 
1ERR*760 , FORM* ' UNFORMATTED ' , READONLY) 

NRECD*1 
GO TO 10 

760 TYPE *, 'CANNOT OPEN FILE* 

GO TO 10 

770. CLOSE (UNIT*8, ERR*780) 

TYPE *, 'CLOSE FILE' 

GO TO 10 

730 TYPE *, 'CANNOT CLOSE FILE' 

GO TO 10 

802 TYPE *, 'FIRST AND LAST RANGE' 

• • ACCEPT * , MRANGE, N RANGE** 

GO TO 10 

803 TYPE 908 / 

ACCEPT * , NITEM yTO PLOT 

GO TO 10 ' 

804 TYPE 911 

ACCEPT *, IHOUR, MINUTE l LEFT EDGE INFO 
FLEFT=3600 . D0* IHOUR+60 . D0*MINUTE 
TYPE 912 

ACCEPT * , XAXIS IFS IN MINUTES 

TYPE 913 

ACCEPT * , IRECDS l# RECORDS TO PLOT 

GO TO 10 

805 TYPE 910 

ACCEPT *,FMULT, SHIFT l PLOT FACTOR 

GO TO 10 

306 TYPE 909 

ACCEPT *» ZERO, FULLS l PLOT SCALE 
GO TO 10 

807 TYPE * , 1 TAPE 13? (0,1) ’ 

ACCEPT * , ITAP13 


808 " 


901 


902 

903 

904 

905 

906 

907 

908 

909 

910 

911 

912 

913 

914 
91.5 

916 

917 

918 

919 

920 

921 

922 
999 


‘3 ) 


GO TO 10 
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TYPE V' COEFFICIENTS 1-4 ( .32, .56, .12, -1.0) : ' 

ACCEPT * , COEF1 , COEF2 , COEF3 , COEF4 
GO TO 10 

FORMAT (IX, 'Xa' ,16, *0 Y= ! ,16, ‘0 N' , 12, ' : • , 12, ' W , 

113, ' s M2, ' FLT ',12, ' RUN ' ,13, ' TIME M2, 

1' f ' , 12, ' : ' ,F4. 1, ' ALT M6) 

FORMAT (IX, 'LOS®' , F6 . 1 , ' THDG® ' , F5 . 1 , ' DRIFT®' , F5 . 2 , 

1' LOSE® ' , F5 . 1 , ’ PITCH® 1 , F6 . 2 , ' ROLL=',F6.2) 

FORMAT (IX, 'TA$i®' , F5 . 1 , * GSPD® ' , F5 . 1 , ' WIND®’,F5.1, 

1' WDIR®' ,F5.i, ' PROC® ' , 06 , ' STATUS®' ,06) 

FORMAT (1X,I5,F6»1,F6.1,I3,I6, '0' , 16 , ' 0 ' , F6 . 2 , ' ',6F6.2) 
FORMAT ( !X, ' AT RECORD', 15,' COMMAND: ',$) 

FORMAT (IX, 'SPACE FORWARD N RECORDS: ' ,$) 

FORMAT (IX, 'SPACE BACKWARD N RECORDS: ',$) 

FORMAT (IX, 'PLOT ITEM M: ',$) 

FORMAT (IX, 'ZERO AND FULL-SCALE FOR PLOT: ',$) 

FORMAT (IX, 'MULTIPLIER AND OFFSET: ' ,$) 

FORMAT (IX, 'LEFT EDGE HOUR, MINUTE: ',$) 

FORMAT (IX, 'FULL SCALE IN MINUTES: ',?) . 

FORMAT (IX, 'NUMBER OF RECORDS: ' ,?) 

FORMAT (IX, 'PLOT CODE TO MATCH (0,8,24)*. ' ,$) ' 

FORMAT (IX, Al, '.(',$) 

FORMAT (IX, Al , $ ) 

FORMAT (IX, 'SPACE TO RECORD N: ',$) 

FORMAT (IX, 5A1 , ' } ' ) 

FORMAT (IX, ' p} ' ,A1, ' . ) ' ) . . . 

FORMAT ( IX, 1 4, 17 ,15 , ' : ' , 12 , ' * ' , F4 . 1, 4F10 .4 ) • •* • ••■ •• • 

FORMAT (IX, 'LO: ' ,F6.2, ' OFFSET: ', F6« 1, * GS:', 3F6.1, 

1' DRIFT: ' , F5 • 2 ) 

FORMAT (IX, 'LOSA®' ,F6.1, ' POWER® ', F6 . 1 , ' VAVG® ' , F6 . 2 , 

1' VCOR® ' , F6 . 2 , ' RECORD® ',15) 

CALL CLOSEC ( ICHAN ) 

STOP 'EXIT 99' 

END 



PROGRAM FED IT 


C 

C 


101 

102 

104 

105 


110 


200 


210 


220 


ORIGINAL PAGE IS 
OF POOR QUALITY 


I NT EGER* 2 HI (56 ) , H2 ( 56 ) , H3 ( 56) , H4( 56) , H5( 56) , 

1D1 ( 1 2 , 90 ) , D2 ( 1 2 , 90 ) , D3 ( 1 2 , 90 ) , D4 (12, 9GT) , D5 (12, 90 ) 
BYTE INPFIL(U),OUTFIL(ll) 

K-0 1 RING BUFFER INDEX 

TYPE 101 

FORMAT (IX, 1 INPUT FILE DL0: * ,$) 

ACCEPT 102, (INPFIL(I) , 1=1,10) 

FORMAT (10A1) 

INPFIL ( 1 1 ) =0 

OPEN ( UNIT-8 , NAME-INPFIL , RECORDS I ZE-568, TYPE- * OLD ' 
1 ERR-900 , FORM- 1 UNFORMATTED ' , READONLY) 

TYPE 105 

FORMAT (IX, 'NUMBER OF RECORDS : ',$) 

ACCEPT * , NRECDS 
IF ( NRECDS. LT.0) GOTO" 104 
NBLK-4 . 47 3*NRECDS+1 
TYPE * , ' NBLK— ' , NBLK 
TYPE 110 

FORMAT (IX, 'OUTPUT FILE DL0:',$) 

ACCEPT 102 , (OUTFIL( I ) , 1=1 , 10 ) 

OUTFIL ( 1 1 ) =0 

OPEN ( UNIT-9 , NAME-OUTFIL , RECORDS! ZE-5.63 , TYPE- ' NEW 
1 ERR-9 10, FORM- ' UNFORMATTED' , I N.ITI ALS I Z E-NBLK )• 

READ ( 8, ERR-920, END-800) H1,D1 
CALL FEDIT1 ( D1 , HI ( 36 ) ) 

READ (8, ERR-920, END-800) H2,D2 
CALL FEDIT1(D2,H2(36) ) 

READ (8, ERR-920, END-800) H3,D3 
CALL FEDIT1(D3,H3(36) ) 

READ (8, ERR-920, END-800) H4,,D4 
CALL FEDIT1(D4,H4(36) ) 

K-K+l 

IF (K.GT.5) K— 1 

GOTO (210,220,230,240,250) K 

STOP 'ILLEGAL K' 

READ (8, ERR-920, END-800) H5,D5 . i 

CALL FEDIT1 ( DS, H5 (36 ) ) 

CALL FEDIT2 (Dl, D3, D5 ) - 
WRITE (9, ERR— 940 ) H3 , D3 
GOTO 200 

READ (8, ERR-920, END-800) H1,D1 
CALL FEDIT1 ( Dl, HI ( 36 ) ) 

CALL FEDIT2 ( D2 , D4 , Dl ) 

WRITE (9,ERR=940) H4,D4 


EDIT FILES CREATED BY FCNVT 
REVISED 10-30-31 AND 11-5-81 
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240 READ (8,ERR=920,END=800) H3,D3 
CALL FEDITi(D3,H3(36) ) 

CALL FEDIT2(D4,D1,D3) 

WRITE (9,ERR=940) H1,D1 
GOTO 200 

250 READ (3,ERR=920,END=800) H4,D4 
CALL FEDIT1 (D4, H4(36 ) ) 

CALL FEDIT2 ( D5, D2 , D4) 

WRITE (9,ERR=940) 112, D2 
GOTO 200 

800 CLOSE (UNIT=8,ERR=930) 

CLOSE (UNIT=9,ERR=930) 

STOP 'COMPLETED* 

900 TYPE *, 'ERROR OPENING INPUT FILE' 
GO TO 800 

910 TYPE *, 'ERROR OPENING OUTPUT FILE' 
GO TO 800 ; 

920 TYPE * , ' READ ERROR* 

GO TO 800 

930 STOP 'CLOSURE ERROR' 

940 TYPE *, 'WRITE ERROR' 

GO TO 800 

END 


230 READ ( 8 , ERR!= s 9^0 , END=800 ) H2,D2 
CALL FEDIT1(D2,H2(36) ) 

CALL FEDIT2(D3,D5,D2) 

WRITE (9,ERR=940) H5,D5 
GOTO 200 


102 


c 

c 


SUBROUT I ME FED XT l ( D, MSNR) 
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LOAD ESTIMATES OF SIGMA INTO D 
REVISED 11-01-81 AND 11-5-81 

INTEGER* 2 D(12, 90) , SIGMA(60) ,MSNR, SNR 

DATA SIGMA/3157, 1571, 1039, 771, 609, 499, 420, 360,313, 
1275,244,218,196,177,161,147,134,123,113,105,97,90, 
184,78,73,68,64,60,57,54,51,48,45,43,41,39,37,35,33, 
132,31,29,28,27,26,25,24,23,22,21,20,20,19,18,18,17, 
116,16,15,15/ 11-SIGMA VEL ERROR PER .2 DB SNR 


DO 200 1=*' 1 , 90 II IS RANGE INDEX 

SNR= ( D ( 1 , I ) -MSNR+200 ) / 20 
IF (SNR.LT.l) SNR=i 
IF (SNR.GT.60) SNR=60 
D(6, I)-SIGMA(SNR) 

200 CONTINUE 


RETURN 

END 
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SUBROUTINE FEDIT2 ( El , E2 ,E3) OF POOR QUALITY 

C PERFORM MAIN EDITING FUNCTION 

C REVISED 11-01-81 AND 11-05-81 

INTEGER* 2 El (12, 90) , E2 (12 ,90) , E3 ( 12, 90) , V(9) , SIGMA, 
1VEL, S(9) ,V1(9) ,VMEAN 

DO 200 IRANGE»1,89 II IS RANGE INDEX 

VEL»E2 ( 2 , I RANGE ) IV, S OF CENTRAL ELEMENT 
SIGMA»E2(6,IRANGE) l CHECK FOR SMALL SIGMA (NO EDIT) 
IF (SIGMA.LT, 40) GOTO 200 

IF (IRANGE.GT.l) GO TO 50 

DO 20 .J-1,3 
JP3=*J+3 
JP6-J+6 
VI ( J)=*E1 ( 2 , J) 

S(J)=E1(6,J) 

V1(JP3)=E2(2, J) 

S(JP3)=E2(6,J) 

V1(JP6)**E3(2,J) 

S(JP6)*E3(6,J) 

20 CONTINUE 
GO TO 100 

50 IM1-IRANGE-1 
. VI(l)=i!l(2,IMl) 

S( 1 )=EiL (6, IM1 ) 

VI ( 2 ) -E2 ( 2 , IM1 ) 

S(2)=E2(6, IMI) 

VI ( 3 ) =E3 ( 2 , IM1 ) 

S( 3 )=E3 (6, IMI) 

IM-IM1+1 
V1(4)=E1(2, IM) 

S(4)=E1 (6, IM) 

VI ( 5 )=E2 ( 2 , IM) 

S(5)=E2 (6, IM) 

VI (6)=E3 ( 2 , IM) 

S(6)=E3(6, IM) 

IP1=IM+1 
VI ( 7 )=E1 ( 2, IP1 ) 

S(7 )=E1 ( 6, IP1) 

VI ( 8 )=E2 ( 2 , IP1 ) 

S(8)=E2(6, IP1) 

VI ( 9)=E3 (2, IP1 ) 

S(9)=E3(6, IP1) 

100 ITHR=300 13 M/S THRESHOLD FOR USE 

105 K=0 
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1 1ST RANGE IS SPECIAL CASE 


! 



»■ 


! 


; jasaMWMg^^ -t£r, 


110 


120 


140 

150 


160 

170 

130 


200 




DO 110 J“l,9 

IF (S(J) .GT . ITHR) GOTO 110 ORIGINAL PAGE IS 
K=K+1 OF POOR QUALITY 

V(K)-V1(J) 

CONTINUE 

IF (K.GT.0) GOTO 120 
ITHR»ITHR*2 

IF ( ITHR. GT. 1200) GOTO 200 
GOTO 105 l TRY AGAIN IF K=0 


M=1 1 SORT V'S THAT PASSED TEST 

IF ( K . LT • 5 ) GOTO 170 INO SORT' REQUIRED 

KM1=K-1 

DO 150 J1=1,KM1 
DO 140 J=1,KM1 

IF (V(J) ,LT.V(J+1)) GOTO 140 
ITEMP»V(J) 

V(J)=V(J+1) 

V(J+1)=ITEMP 

CONTINUE l CAN CONDENSE THIS ROUTINE 

CONTINUE 

MIN=10000 IMINIMIZE V( J)-V( J+3) 

KM3=K*-3 

DO 160 J=1,KM3 
ISUM~V(J)-V(J+3) 

IF (ISUM.GT.MIN) GOTO 160 
MIN-ISUM 

M=J 

CONTINUE.. . • 

IF (K.GT.4) K=4 l FIND MEAN OF MOST LIKELY V’S 

SUM=0 . 

MPKM1=M+K-1 
DO 130 J=M» MPKM1 
SUM=SUM+ FLOAT ( V ( J ) ) 

CONTINUE 

VMEAN=SUM/ FLOAT ( K ) 

I DIFF=I ABS ( VMEAN-VEL ) 

IF (IDIFF.LT. 400) GOTO 200 
IF ( IDIFF. LT . 2*SIGMA) GOTO 200 

E2 ( 7 # IRANGE)=E2 (6 , 1 RANGE) ISAVE FIRST ESTIMATE 

E2(6,IRANGE)=IDIFF-200 

CONTINUE 

RETURN 
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END 
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G SMOOTH FILES CREATED BY FBDIT 

C REVISED 11-2-81 AND 11-5-81 

INTEGER*2 HI ( 56 ) , H2 ( 56 ) , H3 ( 56 ) , H4 ( 56 ) , H5 ( 56 ) , H6 ( 56 ) , H7 ( 56 ) , 
1H8(56),H9(56),D1(12,90) ,J3«I( 12, 90), D3 ( 12,90 ),D4( 12,90 ) , 

1D5 ( 1 2 , 90 ) , D6 ( 1 2 * 90 ) , D7 ( 15,90 ) , D8 ( 1 2 , 90 ) , D9 ( 1 2 , 90 ) 

BYTE INPFIL ( 11 ) , OUTFIL (1 1 ) 

K=0 l RING BUFFER INDEX 

K1=0 

TYPE 101 

101 FORMAT ( IX, ' INPUT FILE DL0:',$) 

ACCEPT 102, (INPFIL(I), 1=1,10) 

102 FORMAT (10A1) 

INPFIL (11) =0 

OPEN ( UNIT=9, NAME=INPFIL, RECORDS I ZE=5 68 , TYPE* ' OLD * , 

1 ERR=900 , FORM= ' UNFORMATTED ' , READONLY) 

104 TYPE 105 

105 FORMAT (IX, 'NUMBER OF RECORDS: ',$) 

ACCEPT * , NRECDS 

IF ( NRECDS. LT.0) GO TO 104 
NBLK=4.473*NRECDS+1 
TYPE * , ' NBLK=' , NBLK 
TYPE 110 

110 FORMAT (IX, 'OUTPUT FILE DL0 : ' , ? ) 

ACCEPT 102, (OUTFIL(I), 1=1,10) 

.. outfil(ii)?=0 :• . .. 

OPEN (UNIT=9, NAME=OUTFIL, RECORDSIZE=568 , TYPE= ' NEW' , 
1ERR=910,FORM=' UNFORMATTED' , INITIALSIZE=NBLK) 

READ (8,ERR=920,END=800) H1,D1 
READ (8,ERR=920,END=800) H2,D2 
READ (8,ERR=920,END=800) H3,D3 
READ (8,ERR=920,END=800) H4,D4 
READ (8, ERR=920 , END=800) H5,D5 
READ ( 8 , ERR=920 , END=800 ) H6 , D6 
READ ( 8 , ERR=920 , END=800 ) H7,D7 
READ (8,ERR=920,END=800) H8,D8 

200 K=K+1 

K1=K1+1 

TYPE * , ' SCAN IN PROCESS :',K1 
IF (K.GT.9) K=1 

GO TO (210,220,230,240,250,260,270,280,290) K 
% STOP 'ILLEGAL K' 

210 READ ( 8 , ERF=920 , END=800 ) H9,D9 
CALL FSMOT1 ( D1 , D3 , D5 , D7 , D9 ) 

WRITE (9,ERR=940) H5,D5 
GO TO 200 

220 READ (8, ERR=920,END=800) H1,D1 
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240 READ (8,ERR=920,END=800) H3,D3 
CALL FSM0T1(D4,D6,D8,D1,D3) 
WRITE ( 9 , ERR«940 ) H8,D8 
GO TO 200 

250 READ ( 8 , ERR**920 , END-800 ) H4,D4 
CALL FSMOTl(D5,D7,D9,D2,D4) 
WRITE (9,ERR*940) H9,D9 
GO TO 200 

260 READ (8,ERR=920,END=800) H5,D5 
CALL FSMOT1 (D6» D8, D,l» D3# D5 ) 
WRITE (9,ERR*940) H1,D1 
GO TO 200 

270 READ (8,ERR=920,EN0s*800) H6,D6 
CALL FSMOT1 (D7 , D9 » D2 , D4# D6 ) 
Write (9,err=940) H2,D2 
GO TO 200 

280 READ ( 8, ERR=920, END=800) H7 , D7 
CALL FSM0T1(D8,D1,D3,D5,D7) 
WRITE ( 9 , ERR=940) H3,D3 
GO TO 200 

» 

290 READ (8, ERR=920 , END=800 ) H8,D8 
CALL FSM0T1(D9,D2,D4,D6,D8) 
WRITE ( 9 , ERR-940 ) H4,D4 
GO TO 200 

800 CLOSE (UNIT=8,ERR=930) 

CLOSE (UNIT=9,ERR=930) 

STOP 'COMPLETED' # 


CALL FSMOTl(D2,D4,D6 f D3,Dl) 
WRITE ( 9 1 ERR*940 ) H6 , D6 
GO TO 200 

230 READ (8, ERR=920, END-800) H2,D2 
CALL FSM0T1(D3,D5,D7,D9,D2) 
WRITE (9, ERR-940) H7,D7 
GO TO 200 


900 

TYPE * , 'ERROR 
GO TO 800 

OPENING INPUT FILE' 

910 

TYPE *, 'ERROR 
GO TO 800 

OPENING OUTPUT FILE' 

920 

TYPE *7 'READ ERROR* 


>G0 TO 800 

. 

930 

STOP ' CLOSURE 

ERROR’ 

940 

TYPE *, 'WRITE 

ERROR' 


GOTO 800 
END 
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SU BROUTINE PS MOT 1(01,02 , 03 , 04 , 05 ) 
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PERFORMS SMOOTHING ON ONE SCAN FROM FSMOTH 
REVISED 03-NO.?'?8l AND 05-NOV-81 

INTEGERS 01 ( 12, 90), D2 (12, 90), 03 (12, 90) , D4( 12, 90) ,05(12,90’/ 
REALM B(25) ,B1(25) ,B2(5) ,6(25, 5) ,G1 (25 , 5 ) ,G2(25 ) ,G4( 15 ) , 
1C1(5) , R(25 ) , W(25 ) , X0, Y0, DX, DY, SUM, SI , S2,T1,T2,T3, STDEV 


INTEGER* 2 IN (5, 5) t POINTER TO SYMMETRIC MATRIX G4 
DATA IN/ l, 2, 4, 7, 11, 2, 3, 5, 8, 12, 4, 5, 6, 
19,13,7,8,9,10,14,11,12,13,14,15/ 

Tl=10000. I THRESHOLDS , 

T2*2 . 5 
T3-1.0 
IER=0 

DO 800 IRANGE=1,87 

IF (IRANGE.GT.3) GO TO 40 
Xii* (IRANGE.EQ.3) GO TO 30 
11*1 
12*5 

IF ( I RANGE. EQ. 2) GO TO 20 

» 

10 ITEST*11 IRANGE*1, LOAD SPECIAL G2 

DATA G2/ . 7 , . 6, . 5 , • 3 , • 2 , » 9 , * 8 , • 8 , » 5 , * 4 , • 9 5 , • 9 , 

1*7, *6, • 5 , .9, .8, *6, . 5 , *4, • 7 , *6, • 5 , • 3 , • 2 / 

GO TO 100 .. ' • 

20 ITEST*12 1 RANGE* 2 , LOAD SPECIAL G2 

DATA G2/ > 6 r • 7 , *6, • 5 , • 3 , • 8 , • 9 , *8, • 8 , *5, * 9 , *95, 

1 • 9 , . 7 , • 6 ^ • 8 , *9, • 8 1 * 6 , *5, • 6 , ♦ 7 , * 8 , *5, < 3/ 

GO TO 100 

30 ITEST=13 1RANGE=3, LOAD NORMAL G2 

DATA G2/ . 5 , • 6 , • 7 , • 6 , *5, * 6 , * 8 , *9 , *8, • 8 , • 7 , *9, 

1*95, . 9, *7, * 6, *8, . 9 , • 3 , *8, *5, *8, *7, *6, *5/ 

40 I 1=1 RANGE-2 l SET UP RANGE OF INDEX 

1 2 *1 RANG E+ 2 • 1 FOR NORMAL RANGES 


100 X0=D3 (4,1 RANGE ) 

Y0=D3 ( 5 , IRANGE ) 
% 

DO 200 1=11,12 


IX, Y OF RELEVANT ELEMENT 

:/ 

l LOAD G, 8 AND W MATRICES 


J1=I-I1+1 l SET UP INDICES 

J2=Jl+5 

J3=J2+5 

J4=J3+5 

J5=J4+5 


B(J1)=D1(2, I)/100. IB IS RAW VELOCITIES 


* 
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200 


B(J2)»D2(2, I)/100. 
B(J3)=D3(2, I)/l00. 
B(J4)»D4(2» I)/l00. 

B ( J 5 ) *D5 (.2,1) /100 * 

W( J1 ) a l00 ./Dl(6 , I) 
W(J2)= S 100./D2(6,I) 

W ( J3 ) » 100 . / D3 ( 6 , 1 ) 
W(J4)»100./D4(6,I) 
W( J5)*100»/D5(6, X) 

DX=*D1(4, l)-X0 
DY»D1 ( 5 , I ) -Y0 
G(J1, 1)=DX**2 
G(Ji,2)=DX 
G ( J1 , 3 ) -DY**2 
G(J1,4)=DY 
G(J1,5)=1. 

DX=*D2(4, I)-X0 
DY=D2(5, 1 ) -Y0 

G(J2,1)-DX**2 
G( J2, 2)=DX 
G(J2,3)=DY**2 
G(J2,4)=DY 
G(J2,5)=l f 

DX=D3(4, I)-X0 
DY=D3(5, l)-Y0 
$(J3,1)=DX**2 
G(J3,2)--DX 
G( J3, 3)= S DY**2 
G(J3,4)=DY 
G ( J3 , 5 ) = 1 . 

DX=D4(4, I)-X0 

DY=D4( 5 , I ) -Y0 

G(J4,1)=DX**2 

G(J4,2)=DX 

G(J4,3)=DY**2 

G(J4,4)=DY 

G ( J4# 5 ) = 1 . 

DX=D5(4, l)-"X0 

DY=D5 (5, I ) -Y0 

G(J5,1)=DX**2 

G(J5,2)=DX 

G(J5,3)=DY**2 

G(J5,4)=DY 

G(J5,5)=1. 

CONTINUE 

IF (W( ITEST) .GT.T3) 
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IW IS 1/ STD DEVtIS 


I FIRST POINT 


l SECOND POINT 


l THIRD POINT 


l FOURTH POINT 


l FIFTH POINT 


W ( ITEST ) =W ( ITEST ) *T2 
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213 SUM*0. I SUM WEIGHTS 

DO 220 1*1,25 
220 SUM=SUM+W(I) 
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IP (SUM.LT.TX) GO TO 300 l EXIT IF THRESHOLD NOT REACHED 


DO 230 1*1,25 

230 W(I)»W(I)*G2(I) l CONVOLVE WEIGHTS 

GO TO 210 

300 CONTINUE 


DO 310 1*1,25 1 WEIGHT LSQ PROBLEM WITH W 

Bl(I)=B(l)*W(l) 

DO 310 J*1 , 5 

310 G 1 ( I , J ) =G (I,J)*W(I) 

DO 320 1*1,5 1 REPLACES GTPRD 

SUM»0. 

DO 315 J=l, 25 

315 SUM=SUM+G 1(J,I)*B1(J) 

320 B2(I\=SUM 

L=0 1G4 IS GlT*Gl 

DO 330 J*l, 5 l THIS ROUTINE REPLACES MATA 

DO 330 1*1, J 
L=L+1 
SUM=0 . 

DO 325 K=1 , 25 

325 SUM=SUM+G.1(R,I)*G1(K,J) . . . 

330 G4(L)=SUM 

CALL SINV(G4,5, .0001, IER) 1G4 IS INVERTED 

DO 360 1*1,5 1 REPLACES MPRD 

SUM=0 *: 

DO 355 J*1 , 5 

355 SUM=SUM+G4 ( IN ( I , J ) ) * B2 ( J ) 

360 Cl ( I)=SUM 


DO 550 1=1,5 

J=I+6 l OUTPUT COEFFICIENTS 

COEF=Cl ( I ) *10000 . 

IF (I.EQ.5) COEF=COEF/100. 1VEL * 1.0 
IP ( COEF.GT. 32767. ) COEF=32767. 

IF ( COEF.LT. -32767. ) COEF=-32767. 

550 D3 ( J, IRANGE)=COEF 

SUM*0 . I FIND MEAN-SQ ERROR IN SOLUTION 

DO 610 1*1,25 

SUM1=-B1(I) 

DO 605 J=l, 5 

605 SUM1=SUM1+G1 ( I, J) *C1 ( J) l G1*C1=RESULTANT VECTOR R * 

610 SUM=SUM+SUMI V * *2 


iio- 


§ 


slBS0> 1 FIND COV MATRIX DIAG SUM 

S2=0. 

DO 620 1=1# 4 

62a S3SwSiStli» tI+l iiiUi« VARIATION • 

S1=S1+ABS(G4(15) ) 

STDEV=SQRT( .2*S1*SUM'I*20>*S2*S2) 
STDEV=STDEV*100. 

IF ( STDEV.GT. 32767 . ) STDEV=32767 . 
D3(12,IRANGE)=STDEV 


800 


CONTINUE 

RETURN 

END 
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PROGRAM FGRID 
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■'£ REDUCE FILES FROM t’SMOTH TO GRIDPOlNTS 
C REVISED NOV 9-1 I 81 

INTEGER*2 HI ( 56 ) , H2 ( 56 ) , H3 ( 56) , H4 ( 56 ) , H5 ( 56 ) , 
1D1(12,90),D2(12,90) ,D3( 12, 90), D4( 12, 90), D5(12,90) 
INTEGER* 2 OUT ( 30 , 20 , 10 ), SOURCE, PAVG 
REAL*4 Z ( 2, 4) 

BYTE INPFIL( 11 ) , OUTFIL( 11 ) 


COMMON OUT, X0, Y0, UCON, VCON, Z, RES, I FLAG, I FORD, 

1 SOURCE , PAVG , ALOS , SCALE , ITHRSH , M RANGE , NRANGE 

NRECD=0 

50 TYPE 52 , NRECD l COMMAND DECODER 

52 FORMAT (IX, 'AT RECORD’, 16,' COMMAND: ',$) 
ACCEPT * , ICOMD 
IF ( ICOMD.LT. 1) GOTO 50 

IF ( IC0MD.GT.7) GOTO 50 I ILLEGAL COMMANDS 

GOTO (100,300,500,600,4*00,800,700) ICOMD 
STOP 'ILLEGAL COMMAND' 


100 

150 

152 

153 

154 

156 

158 

159 

162 


K=0 1 RING BUFFER INDEX 

K1=0 1 EXECUTE GRIDDING OPERATION 

TYPE 150 * 

FORMAT ( IX, ' NUMBER OF SCANS : ' , $ ) 

ACCEPT * , NUMSCN 
TYPE 152 

FORMAT( IX, 'STARTING. X0,Y0 - METERS: ',$)• 
ACCEPT * , X0 , Y0 

X0=X0/20. t **COMMAND LIST** 

Y0=Y0/20. 

TYPE 153 


FORMAT (IX, 'INITIAL AND FINAL RANGE (1-30): ' ,$) 
ACCEPT *,MRANGE, NRANGE 

TYPE 154 11- EXECUTE GRIDDING 


FORMAT (IX, 'AZIMUTH, DEG: '»$) 

ACCEPT *, THETA 
UCON=SIN( THETA/57. 29 58) 
VCON=COS ( THETA/ 57 . 2958 ) 

TYPE 156 1 6 

FORMAT (IX, 'RESOLUTION, METERS 
ACCEPT *, RES 
RES=RES/20. 

TYPE 158 

FORMAT (IX, 'PLOT? (0,1): ',$) 
ACCEPT * , I FLAG 
TYPE 159 


1 2 - OPEN INPUT FILE 
l 3 - SPACE TO RECORD N 
l 4 - DRAW BOX 
i 5 - CLOSE INPUT FILE 
WRITE OUTPUT FILE 
' , $ ) l 7 - STOP PROGRAM 


10 => NO PLOT 


FORMAT (IX, 'SMOOTHED DATA? (0 OR 1):',$) 11=> SMOOTHED DATA 
ACCEPT * , SOURCE 

IF ( I FLAG . LT . 1 ) GOTO 170 l PLOTTER NEEDED? 


TYPE 162 1 SET UP PLOTTER PARMS 

FORMAT (IX, 'M/S PER INCH OF VECTOR: ',$) 
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164 


168 


ACCEPT *, SCALE 
TYPE 164 

FORMAT (IX, 'ERROR THRESHOLD (M/S) s ',$) 
ACCEPT * , THRESH 
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ITHRSH=THRESH*100, 

TYPE 168,155 

FORMAT ( IX, Al, •.(' ,$) 1TURN ON PLOTTER 


170 READ ( 8 , ERR=920 , END=800 ) H1, D1 
READ (9,ERR=920,END=800) H2,D2 
READ ( 8 , ERR=920 , END=800 ) H3,D3 
READ (8,ERR=920,END=800) H4,D4 
NRECD=NRECD+4 


DO 180 11=1,30 
DO 180 12=1,20 
DO 180 13=1,10 
130 OUT (II, 12, 13) =0 


l ZERO OUTPUT MATRIX 


200 K=K+1 l MAIN LOOP RETURN 

IF (K.GT.5) K=1 
K1=K1+1 

IF (K1 .GT.NUMSCN) GOTO 800 
GOTO (210,220,230,240,250) K 
STOP 'ILLEGAL K' 


210 READ ( 8 , ERR=920 , END=800 ) H5,D5 
NRECD=NRECD+1 

IFORD=H3 ( 18 ) I STATUS WORD 

PAVG=H3 ( 36 ) • l MEAN POWER 

ALOS=H3 ( 46 ) 1LOS ANGLE 

CALL FGRID1 (Dl, D3, D5 ) 

GOTO 200 

220 READ ( 8 , ERR=920 , END=800 ) Hl,Dl 
NRECD=NRECD+1 
IFORD=H4( 18 ) 

PAVG=H4 ( 36 ) 

ALOS=H4 ( 46 ) 

CALL FGRID1 ( D2 , D4, Dl ) 

GOTO 200 


230 READ (8, ERR=920, END=800) H2,D2 
NRECD=NRECD+1 
IFORD=H5 ( 18 ) 

PAVG=H5 ( 36 ) 

ALOS=H5 (46 ) 

CALL FGRID1 ( D3 , D5 , D2 ) 

GOTO 200 

240 READ (8, ERR=920,END=800) H3,D3 
NRECD=NRECD+1 
IFORD=Hl( 18 ) 

PAVG=H1 (36 ) 

ALOS=Hl ( 46 ) 
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GOTO 200 !( 

250 READ ( 8 , ERR=920 , END=*800 ) H4,D4 
NRECD=*NRECD+1 
IFORD=H2(18) 

PAVG=H2 ( 36 ) 

ALOS=H2(46) 

CALL FGRID1 ( D5, D2 , D4) 

GOTO 200 

300 TYPE 301 l COMMAND 2 - OPEN FILE 

301 FORMAT (IX, 'INPUT FILE DL0:',$) 

ACCEPT 302 , ( INPFIL( I ) , 1=1, 10) 

302 FORMAT (10A1) 

INPFIL(11)=0 

OPEN ( UNITES , NAME=INPFIL , RECORDSIZE=563 , TYPE® ' OLD ' , 

1 ERR=900 , FORM= ' UNFORMATTED ' , READONLY) 

GOTO 50 

400 CLOSE (UNlT=8,ERR=930) l COMMAND 5 - CLOSE FILE 

,» TYPE *, 'CLOSE INPUT FILE' 

GOTO 50 

500 TYPE 502 V 

502 FORMAT (IX, 'SPACE TO RECORD N: ',$) 

ACCEPT * , NRECD1 

IF ( NRECD1 . EQ . NRECD ) GOTO 50 

IF (NRECD1.LT. 1) GOTO 50 * 

IF (NRECD1.LT. NRECD) GOTO 550 1 BACKSPACE • 
NUMSPC=NRECD1 -NRECD 
DO 520 N=1 , NUMSPC 
READ (8,ERR=920,END=800) H1,D1 
NRECD=NRECD‘f 1 
520 CONTINUE 
GOTO 50 

550 NUMSPC=NRECD-NRECD1 l BACKSPACE REQD 

DO 570 N=l, NUMSPC 
BACKSPACE 8 
NRECD=NRECD-1 
570 CONTINUE 
GOTO 50 

600 TYPE *, 'OUTLINE BOX' 

TYPE 168,155 

CALL GRAPH1(0. ,0. , ' p' ) 

CALL GRAPH1 ( 2000 . , 0 . , ' q ' ) 

CALL GRAPH1( 2000., 1400. ,'q') 

CALL GRAPH1(0., 1400.,' q') 

CALL GRAPH1 ( 0 . , 0 . , 'q* ) . 

TYPE 822,155 
GOTO 50 

700 STOP 'COMPLETED' 
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900 IF (IFLAG.LT. I) GOTO 820 1 PLOTTER ACTIVE? 

TYPE 822,155 

822 FORMAT (IX, ' p) ' ,A1, ' . ) ' ) l TURN OFF PLOTTER 

IFLAG»0 

820 TYPE 801 l WRITE OUTPUT FILE IF REQD 

801 FORMAT ( IX, ' WRITE OUTPUT FILE? (0,1): ’ , $ ) 

ACCEPT * , IWRITE 

IF ( IWRITE. NE.l) GOTO 850 iNO FILE WRITTEN 

TYPE 802 

802 FORMAT ( IX , ' OUTPUT FILENAME DL0 : * , $ ) 

ACCEPT 302, (OUTFIL(I), 1=1,10) 

OUTFIL( 11 )=0 

OPEN ( UNIT=9 , NAME=OUTFIL, RECORDSIZE=3000 , TYPE" ' NEW 
1ERR=910,FORM=' UNFORMATTED' , INITIALSIZE-25 ) 

WRITE (9, ERR=940 ) OUT 
CLOSE ( UNIT=9 , ERR=930 ) 

GOTO 50 

850 TYPE 851 

851 FORMAT (IX, 'PRINT RESULTS (0,1)? ',$) 

ACCEPT *, IWRITE 

IF (I WRITE. NE.l) GOTO 50 INO PRINT RESULTS 
TYPE 855 

855 FORMAT (IX, 'ITEM (1-10): ',$) 

ACCEPT *,NITEM 

TYPE *, 'UI,UF,VI,VF (PRINT LIMITS):' 

ACCEPT *,IUI,IUF,IVI,IVF 

DO 870 IU=IUI , IUF 

PRINT 860, ( OUT ( IU , IV , NITEM ) , IV=IVI , IVF ) 

860 FORMAT (IX, 2016) 

870 CONTINUE 
PRINT * 

PRINT * 

GOTO 50 

900 TYPE *, 'ERROR OPENING INPUT FILE' 

GO TO 50 

910 TYPE *, 'ERROR OPENING OUTPUT FILE' 

GO TO 820 

920 TYPE * , ' READ ERROR' 

GO TO 800 

930 TYPE *, 'CLOSURE ERROR' 

GOTO 50 

940 TYPE * , * WRITE ERROR' 

GO TO 50 


SUBROUTINE FGRID1 ( D1 , D2 , D3 ) 
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REDUCE TO GRIDPOINTS DATA PROM FGRID 
REVISED NOV 9-11 81 

INTEGER*2 D1 ( 12, 90 ) , D2 ( 12 , 90 ) , D3 ( 12 , 90 ) , OUT( 30, 20, 10 ) , 

1 SOURCE, PAVG 
REALM Z ( 2 , 4) 

COMMON OUT ,X0 , Y0 , UCON , VCON , Z , RES , IFLAG, 

1 1 FORD , SOURCE , PAVG , ALOS , SCALE , ITHRSH , M RANGE , N RANGE 
COMMON /LOOP/X, Y 

IFD=0 l FORWARD SCAN? 

IF ( ( IFORD . AND . 24 ) . EQ • 24 ) IFD=5 

DO 800 IRANGE=MRANGE , NRANGE 1 MAIN RANGE LOOP 

IM=IRANGE-1 t TRANSFORM X,Y TO U,V 

IP=IRANGE+1 l AT CORNERS OF 3X3 REGION 

CALL FGRID2(D1(4, IM) ,D1(5, IM) ,1) 

CALL FGRID2 ( D1 (4 , IP j , D1 ( 5 , IP ) , 2 ) 

CALL FGRID2 (D3(4, IM),D3(5,IM),3) 

CALL FGRID2(D3(4, IP) ,D3(5, IP) ,4) 

UMIN=AMIN1 (Z(l,l) ,Z(1,2) ,Z(1,3) ,Z(1,4) ) i FIND EXTREMES OF U, V 
UMAX-AMAX1 ( Z(1 , 1 ) , Z ( 1 , 2 ) , Z ( 1 , 3 ) , Z( 1 , 4) ) 

VMIN=AMIN1 ( Z ( 2 , 1 ) , Z ( 2 , 2 ) , Z ( 2 , 3 ) , Z ( 2 , 4 ) ) 

VM AX=AM AX 1 ( Z( 2 r l ) , Z ( 2 , 2 ) , Z ( 2 , 3.) , Z ( 2 , 4) ) 

IF (UMIN.LT.l.) UMIN=1. i CHECK LEGAL LIMITS 

IF (VMIN.LT.l.) VMIN=1. 

IF (UMAX.LT.l.) GOTO 800 
IF (VMAX.LT.l.) GOTO 800 

IU1=UMIN 1 SET LIMITS OF POSSIBLE GRID PTS 

IV1=VMIN 

IU2=UMAX+1. 

IV2=VMAX+1. 

IF (IU2.GT.30) IU2=30 
IF (IV2.GT.20) IV2=20 

IF (IU1.GT.30) GOTO 800 l IGNORE IF OUT OF' RANGE 

IF (IV1.GT.20) GOTO 800 

***%.■ 

DO 700 IU=IU1 , IU2 IGRIDPOINT LOOP 

DO 700 IV=IV1 , IV2 

X=X0+(IU*UCON-IV*VCON)*RES l GRID POINTS TO X-Y COORDS 
Y=Y0+ ( IU*VCON+I V* UCON ) * RES 

DMIN=FGRID3 ( D2(4 , I RANGE ) , D2 ( 5 , I RANGE) ) 

DO 300 I-IM, IP 
DIST=FGRID3(D1(4, I) , D1 ( 5 , I) ) 

IF (DIST.LT.DMIN) GOTO 700 l SKIP GRIDPOINT IF ANOTHER 


300 

C 


400 


500 


700 

800 


DIST=FGRID3(D2(4, I) , D2(5, I) ) l MEASUREMENT IS CLOSER 
IF (DIST.LT.DMIN) GOTO 700 
DIST=FGRID3 ( D3 (4, I) , D3(5, I) ) 

IF (DIST.LT.DMIN) GOTO 700 

CONTINUE i IF THIS POINT IS REACHED THE POINT AT 

l D2( IRANGE) IS CLOSEST TO GRIDPOINT IU,IV 


OUT(IU, IV, 1+IFD)=D2 ( 1, IRANGE)-PAVG l AMP-NOISE 


0UT(IU, IV, 3+IFD)=D2(3, IRANGE) 
IF ( SOURCE. GT.0) GOTO 400 
OUT ( IU , IV , 2+IFD ) =D2 ( 2 , IRANGE ) 
OUT( IU, IV, 4+IFD)=D2 (6, IRANGE) 


l WIDTH 4 /*—QUT( IU, IV, 5+IFD)*ALO! 
I GOTO 400 IF SMOOTHED 
l RAW VELOCITY 
l EDIT SIGMA 


GOTO 500 


l EXIT TO PLOT TEST 


DX=X-D2 (4, IRANGE) IX, Y DIFFERENCES 

DY=Y-D2 ( 5 , IRANGE ) 

VEL=DX* ( D2( 8, IRANGE )+DX*D2( 7, IRANGE)) l QUADRATIC SURFACE 

VEL=VEL+DY*(D2(10, IRANGE ) +DY*D2 (9, IRANGE) ) 

0UT(IU, IV, 2+IFD) = .0,l*VEL+D2(U, IRANGE) l SMOOTHED VELOCITY 
OUT(IU, IV, 4+IFD)=D2(12, IRANGE) l SMOOTHED SIGMA 


IF (IFLAG.LT.1) GOTO 700 
IF(OUT( IU, IV, 4) . EQ.0 ) GOTO 700 
IF ( OUT ( IU , IV , 9 ) . EQ . 0 ) GOTO 700 
IF(OUT(IU, IV, 4) .GT. ITHRSH) GOTO 700 


IF ( OUT ( IU , IV , 9 ) . GT . ITHRSH ) GOTO 700 


l GOTO 700 IF NO PLOT, 

l OR IF EITHER SIGMA 
IIS ZERO, 

l OR IF AFT OR FORWARD 
I SCANS FAILS TEST 


ALP=>OUT( IU, IV, 5)/572.953 1LOOK ANGLES 

BET=OUT ( IU , IV , 10 ) / 5 7 2 . 958 

Vl=OUT( IU, IV, 2) /100 . I AND RADIAL VELOCITIES 

V2=OUT( IU, IV, 7 ) /100 • 

SINBMA=SIN ( BET-ALP ) 

XX=(Vl*COS(BET)-V2*COS(ALP) )/SINBMA !X,Y FOR PHI=90 DEG 

YY= ( V2*SIN ( ALP) -V1*SIN( BET) ) / SINBMA 
UU=YY*VCON+XX*UCON IU,V COORDS IN M/S 

W=YY* UCON-XX*VCON 

DX=100.*UU/ SCALE 
DY=100.*W/ SCALE 
IF (ABS(DX) .GT.400. ) GOTO 
IF (ABS(DY) .GT.400. ) GOTO 
X=60.*IU+100.-DX 
Y=60.*IV+100.-DY 
CALL GRAPH1 (X, Y, 'p 1 ) 

X=X+2.*DX 
Y=Y+2.*DY 

CALL GRAPH1 (X, Y, ' q' ) 

X=X-.5*DX-.2*DY 
Y=Y- . 5*DY+ . 2*DX 
CALL GRAPH1(X,Y, * q ’ ) 
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CONTINUE 

CONTINUE 

RETURN 

END 


l VECTOR HALF-LENGTHS 

700 l EXIT IF VECTOR TOO LONG 

700 

! ORIGIN OF VECTOR 

l MOVE TO ORIGIN 

l PLOT BODY OF VECTOR 

I PLOT BARB 
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SUBROUTINE FGRID2(X, Y, I) 

TRANSFORM X,Y COORDINATES INTO U,V COORDINATES 
REVISED 9-NOV-81 

INTEGER* 2 OUT( 30 , 20, 10 ) , SOURCE, PAVG,X, Y 
REAL *4 Z ( 2 , 4 ) 

COMMON OUT , X0 , Y0 , UCON , VCON , Z , RES , I FLAG, I FORD, 

1 SOURCE , PAVG , ALOS , SCALE , ITHRSH 

Z (1 , 1 ) - ( ( X-X0 ) * UCON+ ( Y-Y0 ) * VCON ) / RES 
Z(2,I)=( ( Y-Y0 ) *UCON- ( X-X0 ) * VCON ) / RES 
RETURN 

END 


REAL FUNCTION FGRID3 (XI , Yi ) 

CALCULATES SQUARE-DISTANCE TO X,Y 
CREATED 10-NOV-81 

INTEGER*2 XI, Yl 
COMMON /LOOP/X, Y 

FGRID3=(X1-X)**2+(Y1-Y)**2 

RETURN 

END 


It 






*• 


♦ 


it* 


950 


.. • m • BA/iir 


SUBROUTINE GRAPH! (X, Y, A) 

BYTE OUT(5) , A 

IF (X.GT. 32767.) X=32767. 

IF (Y.GT. 32767.) Y=32767. 

IF (X.LT.0.) X=0. 

IF (Y.LT.0. ) Y=0. 

IX1=X/ 1024 

IX2=X/l6-64*IXl 

1X3 =X-1 6 *1X2 -102 4*1X1 

IYl=Y/4096 

IY2=Y/64-64*IYl 

OUT( 1 ) =224+1X1 

OUT( 2) =128+1X2 

IF (0UT(2) .LT.160) OUT ( 2 ) --OUT ( 2 ) +64 
OUT ( 3 ) =128+IY1+4*IX3 
IF (OUT(3) .LT.160) OUT(3 ) =OUT( 3 ) +64 
OUT ( 4 ) =1 28+IY2 

IF (OUT(4) .LT.160) OUT(4)=OUT(4)+64 
OUT (5 )=( 128 .+Y-64 . *IY2-4096 . *IY1 ) 

IF (OUT (5) .LT.160) OUT ( 5 ) =OUT ( 5 ) +64 
TYPE 950, A, OUT 
FORMAT ( IX, 6A1 , ' } ' ) 

RETURN 



END 



t 
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1 


2 


3 

4 

5 

6 


7 

8 
9 
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l * , jr 


SUBROUTINE SINV( A, N, EPS* IER) 
DIMENSION A(I) 

REAL*8 DIN, WORK 
CALL MFSD(A*N* EPS, IER) 

IF(IER) 9*1*1 
IPIV»N*(N+l)/2 
IND-IPIV 
DO 6 1=1, N 

DIN=1.0D0/DBLE(A(IPIV) ) 

A( IPIV) =DIN 

MIN-N 

KEND=I-1 

LANF=N-KEND 

IF(KEND) 5,5,2 

J=IND 

DO 4 K=1 , KEND 
WORK=0 . 0D0 
MIN=MIN-1 
LHOR=IPIV 
LVERwJ 

DO 3 L=LANF, MIN 

LVER=LVER+1 

LHOR=LHOR+L 

WORK=WORK+DBLE ( A(LVER) *A( LHOR) ) 

A( J)=-WORK*DIN 

J=J-MIN 

IPIV=IPIV-MIN 

IND=IND-1 

CONTINUE 

DO 8 1=1, N 

IPIV=IPIV+I 

J=IPIV 

DO 8 K=I , N 

WORK=0 . 0D0 

LHOR=J 

DO 7 L=K, N 

LVER=LHOR+K- J. 

WORK=WORK+DBLE ( A( LHOR) *A( LVER) ) 

LHOR=LHOR+L 

A(J)=WORK 

J=J+K 

RETURN 

END 
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* 


41 ? 


* 
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SUBROUTINE MFSD( A, N, EPS, IER) 
DIMENSION A( 1 ) 

REAL*8 DPIV , DSUM 

IP (N-l ) 12,1,1 

IER-0 

KPIV-0 

DO 11 K=*l , N 

KPIV=KPIV+K 

IND=KPIV 

LEND=K-1 

TOL*ABS(EPS*A(KPIV) ) 

DO 11 I=K,N 
DSUM=0.0D0 
IF (LEND) 2,4,2 
DO 3 L=1 , LEND 
LANF=KPIV-L 
LIND=IND-L 

DSUM=DSUM+DBLE ( A( LANF ) * A ( LIND ) ) 

DSUM-DBLE ( A( IND ) ) -DSUM 

IP(I-K) 10,5,10 

IF ( SNGL ( DSUM) -TOL) 6,6,9 

IF ( DSUM) 12,12,7 

IF ( IER) 8,8,9 

IER=K-1 

DPIV=DSQRT( DSUM) 

A(KPIV)=DPIV 
DPIV=1 . 0D0/DPIV 
GO TO 11 

A( IND) =DSUM*DPIV 

IND=IND+I 

RETURN 

IER=-1 

RETURN 

END 
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